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Abstract 

We construct a family of stochastic growth models in 2 + 1 dimen- 
sions, that belong to the anisotropic KPZ class. Appropriate projec- 
tions of these models yield 1 + 1 dimensional growth models in the 
KPZ class and random tiling models. We show that correlation func- 
tions associated to our models have determinantal structure, and we 
study large time asymptotics for one of the models. 

The main asymptotic results are: (1) The growing surface has 
a limit shape that consists of facets interpolated by a curved piece. 
(2) The one-point fluctuations of the height function in the curved 
part are asymptotically normal with variance of order ln(i) for time 
t 3> 1. (3) There is a map of the (2 + l)-dimensional space-time to the 
upper half-plane EI such that on space-like submanifolds the multi- 
point fluctuations of the height function are asymptotically equal to 
those of the pullback of the Gaussian free (massless) field on HI. 
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1 Introduction 



In recent years there has been a lot of progress in understanding large time 
fluctuations of driven interacting particle systems on the one- dimensional 
lattice, see e.g. [1,2,4,7-11,25,29,30,38,45-48]. Evolution of such systems 
is commonly interpreted as random growth of a one-dimensional interface, 
and if one views the time as an extra variable, the evolution produces a 
random surface (see e.g. Figure 4.5 in [43] for a nice illustration). In a 
different direction, substantial progress have also been achieved in studying 
the asymptotics of random surfaces arising from dimers on planar bipartite 
graphs, see the review [32] and references therein. 

Although random surfaces of these two kinds were shown to share certain 
asymptotic properties (also common to random matrix models), no direct 
connection between them was known. One goal of this paper is to establish 
such a connection. 

We construct a class of two-dimensional random growth models (that is, 
the principal object is a randomly growing surface, embedded in the four- 
dimensional space-time). In two different projections these models yield ran- 
dom surfaces of the two kinds mentioned above (one reduces the spatial 
dimension by one, the second projection is fixing time). We partially com- 
pute the correlation functions of an associated (three-dimensional) random 
point process and show that they have determinantal form that is typical for 
determinant al point processes. 

For one specific growth model we compute the correlation kernel explic- 
itly, and use it to establish Gaussian fluctuations of the growing random 
surface. We then determine the covariance structure. 

Let us describe our results in more detail. 

1.1 A two-dimensional growth model 

Consider a continuous time Markov chain on the state space of interlacing 
variables 

= mCZ^I^^^/^r}. n = l,2,.... (1.1) 

m=l,...,n ' 

x™ can be interpreted as the position of particle with label (k,m), but we 
will also refer to a given particle as x™. As initial condition, we consider the 
fully-packed one, namely at time moment t = we have x™(0) = k — m — 1 
for all k, m, see Figure 11.11 

The particles evolve according to the following dynamics. Each of the 
particles x™ has an independent exponential clock of rate one, and when the 



3 




Figure 1.1: Illustration of the initial conditions for the particles system and 
the corresponding lozenge tilings. In the height function picture, the white 
circle has coordinates (x,n,h) = (—1/2,0,0). 



x™-clock rings the particle attempts to jump to the right by one. If at that 
moment x™ = x™~ x — 1 then the jump is blocked. If that is not the case, 
we find the largest c > 1 such that x™ = = ■ ■ ■ = x™^-i > an d all c 

particles in this string jump to the right by one. For any t > denote by 
Ai^ n \t) the resulting measure on at time moment t. 

Informally speaking, the particles with smaller upper indices are heavier 
than those with larger upper indices, so that the heavier particles block and 
push the lighter ones in order for the interlacing conditions to be preserved. 
This anisotropy is essential, see more details in Section 11.41 

Let us illustrate the dynamics using Figure 11.21 which shows a possible 
configuration of particles obtained from our initial condition. If in this state 
of the system the a^-clock rings, then particle x\ does not move, because it is 
blocked by particle x\. If it is the x|-clock that rings, then particle x\ moves 
to the right by one unit, but to keep the interlacing property satisfied, also 
particles x\ and x\ move by one unit at the same time. This aspect of the 
dynamics is called "pushing". 

Observe that «S (ni) C <S ( ™ 2) for 

n i < n 2, an d the definition of the evolution 
implies that Ai^ ni \t) is a marginal of Ai^ U2 \t) for any t > 0. Thus, we can 
think of .M^'s as marginals of the measure Ai = lim Ai^ on S = limS^ n \ 

In other words, Ai(t) are measures on the space S of infinite point configu- 
rations W}*=l,...,m >m > r 

Before stating the main results, it is interesting to notice that the Markov 
chain has different interpretations. Also, some projections of the Markov 
chain to subsets of are still Markov chains. 
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Charlier process 




Figure 1.2: From particle configurations (left) to 3d visualization via lozenge 
tilings (right). The corner with the white circle has coordinates (x,n,h) = 
(-1/2,0,0). 



f . The evolution of x\ is the one-dimensional Poisson process of rate one. 

2. The row {x™} m >i evolves as a Markov chain on Z known as the Totally 
Asymmetric Simple Exclusion Process (TASEP), and the initial condi- 
tion x™(0) = — m is commonly referred to as step initial condition. In 
this case, particle x\ jumps to its right with unit rate, provided the 
arrival site is empty (exclusion constraint). 

3. The row {:r™} m >i also evolves as a Markov chain on Z that is sometimes 
called "long range TASEP"; it was also called PushASEP in [7]. It is 
convenient to view {x™ + m} m >i as particle locations in Z. Then, 
when the x^-clock rings, the particle x\ + k jumps to its right and 
pushes by one unit the (maybe empty) block of particles sitting next 
to it. If one disregards the particle labeling, one can think of particles 
as independently jumping to the next free site on their right with unit 
rate. 

4. For our initial condition, the evolution of each row {x™}k=i,...,m, 
m — 1,2, ... , is also a Markov chain. It was called Charlier process 
in [37] because of its relation to the classical orthogonal Charlier poly- 
nomials. It can be defined as Doob /i-transform for m independent 
rate one Poisson processes with the harmonic function h equal to the 
Vandermonde determinant. 

5. Infinite point configurations {x™} G S can be viewed as Gelfand-Tsetlin 
schemes. Then M. (t) is the "Fourier transform" of a suitable irreducible 
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character of the infinite-dimensional unitary group U(oo), see [13]. In- 
terestingly enough, increasing t corresponds to a deterministic flow on 
the space of irreducible characters of U(oo). 

6. Elements of S can also be viewed as lozenge tiling of a sector in the 
plane. To see that one surrounds each particle location by a rhombus 
of one type and draws edges through locations where there are no par- 
ticles, see Figure 11.21 Our initial condition corresponds to a perfectly 
regular tiling, see Figure 11.11 

7. The random tiling defined by Ai(t) is the limit of the uniformly dis- 
tributed lozenge tilings of hexagons with side lengths (a, b, c), when 
a, b, c — > oo so that ab/c — > t, and we observe the hexagon tiling at 
finite distances from the corner between sides of lengths a and b. 

8. Finally, Figure 11.21 has a clear three-dimensional connotation. Given 
the random configuration {x^(t)} 6 <S at time moment t, define the 
random height function 

h : (Z + |) x Z >0 x M> -> Z> , 

h(x,n,t)=#{ke{l,...,n}\x%t)>x}. [ ' ] 

In terms of the tiling on Figure IT!2l the height function is defined at the 
vertices of rhombi, and it counts the number of particles to the right 
from a given vertex. (This definition differs by a simple linear function 
of (x, n) from the standard definition of the height function for lozenge 
tilings, see e.g. [32,33].) The initial condition corresponds to starting 
with perfectly flat facets. 

Thus, our Markov chain can be viewed as a random growth model of the 
surface given by the height function. In terms of the step surface of Figure [I~2l 
the evolution consists of removing all columns of (x, n, h) -dimensions (1, *, 1) 
that could be removed, independently with exponential waiting times of rate 
one. For example, if x\ jumps to its right, then three consecutive cubes 
(associated to removed. Clearly, in this dynamics the directions 

x and n do not play symmetric roles. Indeed, this model belongs to the 2 + 1 
anisotropic KPZ class of stochastic growth models, see Section 11.41 
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1.2 Determinantal formula, limit shape and one-point 
fluctuations 

The first result about the Markov chain Ai(t) that we prove is the (partial) 
determinantal structure of the correlation functions. Introduce the notation 

(ni,ti) -< (n 2 ,t 2 ) iff rti < n 2 ,t 1 > t 2 , and (n l ,ti) ^ (n 2 ,t 2 ). (1.3) 

Theorem 1.1. For any N = 1, 2, . . . , pick N triples 

Kj = (Xj,7lj,tj) G Z x Z >0 x R> 

such that 

ti < t 2 < ■ • • < t N , n\ > n 2 > ■ • • > n N . (1.4) 

Then 

P{For each j = 1, . . . , N there exists a kj, 

1 < kj < rij such that x^ifj) = Xj} = det [/C(xj, Xj)]fj =v (1.5) 

where 

1 r dw e (*i-t2)A» 

/C(Xi, n h h] X 2 , n 2 , t 2 ) — w x 2 -x 1+ l ( X _ w )n 2 -m 1 [(ni,ti)-<(n2,ta)] 

1 / , / , (1 - w)™ 1 w Xl 1 . . 

+ 7 777 dw ® dz i , - — , 1.6 

(27ri) 2 J To J Tl e^l z (1 - z) n * w - ;z v ; 

£/ie contours T$, T\ are simple positively oriented closed paths that include 
the poles and 1, respectively, and no other poles (hence, they are disjoint). 

This result is proved at the end of Section |2~%1 The above kernel has in fact 
already appeared in [7] in connection with PushASEP. The determinantal 
structure makes it possible to study the asymptotics. On a macroscopic scale 
(large time limit and hydrodynamic scaling) the model has a limit shape, 
which we now describe, see Figure 11.31 Since we look at heights at different 
times, we cannot use time as a large parameter. Instead, we introduce a large 
parameter L and consider space and time coordinates that are comparable 
to L. The limit shape consists of three facets interpolated by a curved piece. 
To describe it, consider the set 

V = {(u, 77, r) G M 3 >0 I (VV ~ v^) 2 < v < (y/Tj + v^) 2 }. (1.7) 

It is exactly the set of triples (y, rj, r) G R> for which there exists a non- 
degenerate triangle with side lengths (y/v, y/fj, Vt). Denote by (tt u , n v , tc t ) 
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the angles of this triangle that are opposite to the corresponding sides (see 
Figure EH] too). 

Our second result concerns the limit shape and the Gaussian fluctuations 
in the curved region, living on a vhiL scale. 

Theorem 1.2. For any (i/,t),t) G T> we have the moment convergence of 
random variables 

hm Mt^-^ + i^.r^-EMt^-^ + i^rL) = ^ 

(1.8) 

tflitb « = (27T 2 )" 1 . 

We also give an explicit formula for the limit shape: 

hm =: h(u,r],T) 

If , , sin 7i u sin tt„ \ 
= -l-un ri + r l (n-n u ) + T ^ v - . 1.9 

7T \ Sin 7T T J 

Theorem 11.21 describes the limit shape h of our growing surface, and 
the domain T> describes the points where this limit shape is curved. The 
logarithmic fluctuations is essentially a consequence of the local asymptotic 
behavior being governed by the discrete sine kernel (this local behavior occurs 
also in tiling models [24,31,42]). Using the connection with the Charlier 
ensembles, see above, the formula (11.91) for the limit shape can be read off 
the formulas of [5]. 

Using Theorem 11.11 it is not hard to verify (see Proposition 13.21 below) 
that near every point of the limit shape in the curved region, at any fixed 
time moment the random lozenge tiling approaches the unique translation 
invariant measure M Wvt7T %T on lozenge tilings of the plane with prescribed 
slope (see [16,32,35] and references therein for discussions of these measures). 
The slope is exactly the slope of the tangent plane to the limit shape, given 
by 

av 7r or] 7r 

This implies in particular, that (tTu/tt, tt^/it, ti t /ti) are the asymptotic pro- 
portions of lozenges of three different types in the neighborhood of the point 
of the limit shape. One also computes the growth velocity (see flTTTZJ) for the 
definition of Q) 

dh 1 sin tt v sin ir v lm(Q(u, 77, r)) (111) 
dr ix sin 71V vr 
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Since the right-hand side depends only on the slope of the tangent plane, 
this suggest that it should be possible to extend the definition of our surface 
evolution to the random surfaces distributed according to measures M 7r ^ i7r , ;)7rT ; 
these measures have to remain invariant under evolution, and the speed of 
the height growth should be given by the right-hand side of ( 11. lip . This is 
an interesting open problem that we do not address in this paper. 



1.3 Complex structure and multipoint fluctuations 

To describe the correlations of the interface, we first need to introduce a com- 
plex structure. Set EI = {z 6 C | lm(z) > 0} and define the map Q : V — ► EI 
by 

p(u,rj,r)\ = v 7 ^, \l-n(u,r],r)\ = ^Tfr. (1.12) 

Observe that argf2 = tt u and arg(l — Q) — —7i n . The preimage of any Q e H 
is a ray in T> that consists of triples (u, r),r) with constant ratios {y : 77 : r). 
Denote this ray by Rq. One sees that ffo's are also the level sets of the slope 
of the tangent plane to the limit shape. Since h{av, arj, ar) = ahiv^ rj, r) for 
any a > 0, the height function grows linearly in time along each Rq. Note 
also that the map Q satisfies 

and the first of these relations is the complex Burgers equation, cf . [34] . 

From Theorem 11.21 one might think that to get non-trivial correlations 
we need to consider {h — E(/i))/VhiL. However, this is not true and the 
division by VhiL is not needed. To state the precise result, denote by 

z — w 



Q{z,w) = -— In 
In 



w 



1.14) 



the Green function of the Laplace operator on EI with Dirichlet boundary 
conditions. 

Theorem 1.3. For any N = 1, 2, . . . , let Xj = (i/j, rjj, Tj) ET> be any distinct 
N triples such that 

n < T 2 < ■ ■ ■ <T N , 771 > 7] 2 > ■ ■ ■ > Tj N . (1.15) 

Denote 

H L (u, r), r) := v^F (h{[{u - r,)L] + §, [rjL], tL) - Eh([(v - rj)L) + ±, [t]L],tL)) 

(1.16) 
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and flj = Q(vj,T]j, Tj). Then 



lim E{H L {x 1 )---H L {x N )) 

L— >oo 



N/2 

E II ^(^<x(2i-i),fi<7(2j)), N is even, 

o^Tn 3=1 



0, zs odd, 

(1-17) 

where the summation is taken over all fixed point free involutions a on 
{1 V}. 

The result of the theorem means that as L — * oo, i7^(n _1 (z)) is a Gaus- 
sian process with covariance given by Q, i.e., it has correlation of the Gaus- 
sian Free Field on EL We can make this statement more precise. Indeed, 
in addition to Theorem 11.31 a simple consequence of Theorem 11.21 gives (see 
Lemma 15.41) , 

E ■ • ■ H L {x N )) = 0(L e ), L -> oo, (1.18) 

for any Kj G T> and any e > 0. This bounds the moments of Hl(xj) for 
infmitesimally close points Xj. A small extension of Theorem 11.31 together 
with this estimate immediately implies that on suitable surfaces in T>, the 
random function Hl(v, rj, r) converges to the f2-pullback of the Gaussian free 
field on H, see Theorem 15.61 and Theorem 15.81 in Section I5~5l for more details. 

Conjecture 1.4. The statement of Theorem \1.3\ holds without the assump- 
tion U.15\) . provided that Q-images of all the triples are pairwise distinct. 

Theorem 11.31 and Conjecture 11.41 indicate that the fluctuations of the 
height function along the rays Rq vary slower than in any other space-time 
direction. This statement can be rephrased more generally: the height func- 
tion has smaller fluctuations along the curves where the slope of the limit 
shape remains constant. We have been able to find evidence for such a claim 
in one-dimensional random growth models as well [23]. 



1.4 Universality class 

In the terminology of physics literature, see e.g. [3], our Markov chain falls 
into the class of local growth models with relaxation and lateral growth, 
described by the Kardar-Parisi- Zhang (KPZ) equation 

d t h = Ah + Q(d x h, d y h) + white noise, (1-19) 

where Q is a quadratic form. Relations fll.lOp and (11. lip imply that for our 
growth model the determinant of the Hessian of dth, viewed as a function of 
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Figure 1.3: A configuration of the model analyzed with N = 100 particles at 
time t = 25, using the same representation as in Figure ll.2[ In [22] there is 
a Java animation of the model. 



the slope, is strictly negative, which means that the form Q in our case has 
signature (—1, 1). In such a situation the equation (11.191) is called anisotropic 
KPZ or AKPZ equation. 

An example of such system is growth of vicinal surfaces, which are nat- 
urally anisotropic because the tilt direction of the surface is special. Using 
non-rigorous renormalization group analysis based on one-loop expansion, 
Wolf [52] predicted that large time fluctuations (the roughness) of the growth 
models described by AKPZ equation should be similar to those of linear 
models described by the Edwards- Wilkinson equation (heat equation with 
random term) 

d t h = Ah + white noise. (1-20) 

Our results can be viewed as the first rigorous analysis of a non- 
equilibrium growth model in the AKPZ class. (Some results, like loga- 
rithmic fluctuations, for an AKPZ model in a steady state were obtained 
in [44]. Some numerical numerical results are described in [27,28,36]). In- 
deed, Wolf's prediction correctly identifies the logarithmic behavior of height 
fluctuations. However, it does not (at least explicitly) predict the appearance 
of the Gaussian free field, and in particular the complete structure (map Q) 
of the fluctuations described in the previous section. 

On the other hand, universality considerations imply that analogs of The- 
orems 11.21 and 11.31 as well as possibly Conjecture 11.41 should hold in any 
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AKPZ growth model. 

1.5 More general growth models 

It turns out that the determinantal structure of the correlations functions 
stated in Theorem 11.11 holds for a much more general class of two-dimensional 
growth models. In the first part of the paper we develop an algebraic for- 
malism needed to show that. At least three examples where this formalism 
applies, other than the Markov chain considered above, are worth mention- 
ing. 

1. In the Markov chain considered above one can make the particle jump 
rates depend on the upper index m in an arbitrary way. One can also 
allow the particles jump both right and left, with ratio of left and right 
jump rates possibly changing in time [7]. 

2. The shuffling algorithm for domino tilings of Aztec diamonds intro- 
duced in [21] also fits into our formalism. The corresponding discrete 
time Markov chain is described in Section [2] below, and its equivalence 
to domino shuffling is established in the recent paper [39]. 

3. A shuffling algorithm for lozenge tilings of the hexagon (also known as 
boxed plane partitions) has been constructed in [12] using the formalism 
developed in this paper, see [12] for details. 

Our original Markov chain is a suitable degeneration of each of these 
examples. 

We expect our asymptotic methods to be applicable to many other two- 
dimensional growth models produced by the general formalism, and we plan 
to return to this discussion in a later publication. 

1.6 Other connections 

We have so far discussed the global asymptotic behavior of our growing sur- 
face, and its bulk properties (measures M nv>% _ l7rT ), but have not discussed the 
edge asymptotics. As was mentioned above, rows {x™} m >i and {x^}m>i can 
be viewed as one-dimensional growth models on their own, and their asymp- 
totic behavior was studied in [7] using essentially the same Theorem 11.11 
This is exactly the edge behavior of our two-dimensional growth model. 

Of course, the successive projections to {x™} m >i and then to a fixed 
(large) time commute. In the first ordering, this can be seen as the large time 
interface associated to the TASEP. In the second ordering, it corresponds to 
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considering a tiling problem of a large region and focusing on the border of 
the facet. 

Interestingly enough, an analog of Theorem 11.11 remains useful for the 
edge computations even in the cases when the measure on the space S is no 
longer positive (but its projections to {x™}m>i and {x™}m>i remain pos- 
itive). These computations lead to the asymptotic results of [7-11,48] for 
one-dimensional growth models with more general types of initial conditions. 

Another natural asymptotic question that was not discussed is the lim- 
iting behavior of MS n \t) when t — ► oo but n remains fixed. After proper 
normalization, in the limit one obtains the Markov chain investigated in [51]. 

Two of the four one-dimensional growth models constructed in [20] 
(namely, "Bernoulli with blocking" and "Bernoulli with pushing" ) are projec- 
tions to {x™} m >i and {x™}m>i of one of our two-dimensional growth models, 
see Section [2] below. It remains unclear however, how to interpret the other 
two models of [20] in a similar fashion. 

Finally, let us mention that our proof of Theorem 11.11 is based on the ar- 
gument of [18] and [50], the proof of Theorem [0] uses several ideas from [33], 
and the algebraic formalism for two-dimensional growth models employs a 
crucial idea of constructing bivariate Markov chains out of commuting uni- 
variate ones from [19]. 

Outline. The rest of the paper is organized as follows. It has essentially 
two main parts. The first part is Section [2j It contains the construction of 
the Markov chains, with the final result being the determinantal structure 
and the associated kernel (Theorem 12.25ft . Its continuous time analogue is 
Corollary 12.261 whose further specialization to particle-independent jump 
rate leads to Theorem 11.11 The second main part concerns the limit results 
for the continuous time model that we analyze. We start by collecting various 
geometric identities in Section [3J We also shortly discuss why our model is 
in the AKPZ class. In Section H] we first give a shifted version of the kernel, 
whose asymptotic analysis is the content of Section [UJ These results then 
allow us to prove Theorem 11.21 in Section H] and Theorem 11.31 in Section [5j 

Acknowledgments. The authors are very grateful to P. Diaconis, 
E. Rains, and H. Spohn for numerous illuminating discussions. The first 
named author (A. B.) was partially supported by the NSF grant DMS- 
0707163. 
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2 Two dimensional dynamics 



All the constructions below are based on the following basic idea. Consider 
two Markov operators P and P* on state spaces S and S*, and a Markov 
link A : S* — > S that intertwines P and P*, that is AP = P*A. Then one 
can construct Markov chains on (subsets of) S* x S that in some sense has 
both P and P* as their projections. There is more than one way to realize 
this idea, and in this paper we discuss two variants. 

In one of them the image {y*,y) of (x*,x) G S* x S under the Markov 
operator is determined by sequential update: One first chooses y according 
to P(x,y), and then one chooses y* so that the needed projection properties 
are satisfied. A characteristic feature of the construction is that x and y* are 
independent, given x* and y. This bivariate Markov chain is denoted Pa; its 
construction is borrowed from [19]. 

In the second variant, the images y* and y are independent, given (x, x*), 
and we say that they are obtained by parallel update. The distribution of y 
is still P(x ) y) ) independently of what x* is. This Markov chain is denoted 
Pa for the operator A = AP = P*A that plays an important role. 

By induction, one constructs multivariate Markov chains out of finitely 
many univariate ones and links that intertwine them. Again, we use two 
variants of the construction — with sequential and parallel updates. 

The key property that makes these constructions useful is the following: 
If the chains P, P*, and A, are /i-Doob transforms of some (simpler) Markov 
chains, and the harmonic functions h used are consistent, then the transition 
probabilities of the multivariate Markov chains do not depend on h. Thus, 
participating multivariate Markov chains may be fairly complex, while the 
transition probabilities of the univariate Markov chains remain simple. 

Below we first explain the abstract construction of Pa, Pa, and their mul- 
tivariate extensions. Then we exhibit a class of examples that are of interest 
to us. Finally, we show how the knowledge of certain averages (correlation 
functions) for the univariate Markov chains allows one to compute similar 
averages for the multivariate chains. 

2.1 Bivariate Markov chains 

Let S and S* be discrete sets, and let P and P* be stochastic matrices on 
these sets: 

£P(s,y) = l, xeS; Yl P*(x m ,V m ) = l, e (2-1) 
Assume that there exists a third stochastic matrix 
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A = ||A(x* 

' x ) Wx*eS* xes suc h- that for any x* G S* and y G S 

J2Hx*,x)P(x,y) = P*{x\y*)K{y\y). (2.2) 
xes y*es* 

Let us denote the above quantity by A(x*,y). In matrix notation 

A = AP = P*A. (2.3) 

Set 

S A = {( x *, x) G S* x S | A(x*, x) > 0}, 
S A = {{x*, x) G S* x 5 | A(x*, rr) > 0}. 

Define bivariate Markov chains on S\ and Sa by their corresponding transi- 
tion probabilities 

( P(x,y)P*(x*,y*)A(y*,y) 
P A ((x*,x),(y*,y))=l A(x*,y) ~ ' ^ > U ' (2 .4) 

[0, otherwise, 

PM ^ x) , w , y)) = ^v}^pLsL. (2 .5) 

It is immediately verified that both matrices Pa and Pa are stochastic. 

The chain Pa was introduced by Diaconis-Fill in [19], and we are using 
the notation of that paper. 

One could think of Pa and Pa as follows. 

For Pa, starting from (x*,x) we first choose y according to the transition 
matrix P(x,y), and then choose y* using P ^^J^ ' y \ which is the con- 
ditional distribution of the middle point in the successive application of P* 
and A provided that we start at x* and finish at y. 

For Pa, starting from (x*,x) we independently choose y according to 
P(x,y) and y* according to P ^^jj^ ' x \ which is the conditional distribu- 
tion of the middle point in the successive application of P* and A provided 
that we start at x* and finish at x. 

Lemma 2.1. For any (x*,x) G S A , y G S, we have 



P A ((x*,x),(y*,y)) = P(x,y), 

y*eS*:(y*,y)eS A 

P A ((x*,x),(y*,y)) = P(x,y), 

y*eS*:(y*,y)eS A 



(2.6) 
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and for any f £ 5*, (y*, y) G 5a, 

A(x*,x)P A ((x*, x), (y*, y)) = P*(x*, y*)A(y*, y), 
xeS:(x*,x)es A ^2 ^ 

A(x*, x)P A {(x*, x), {y*, y)) = P*(x*, y*) A(y\ y). 

Proof of Lemma \2.1\ Straightforward computation using the relation 
A = AP = P*A. □ 

Proposition 2.2. Let m*(x*) be a probability measure on S* . Consider the 
evolution of the measure m(x*)A(x* ,x) on S\ under the Markov chain P\ 
and denote by (x*(j),x(j)) the result after j = 0,1,2,... steps. Then for 
any k, I = 0, 1, . . . the joint distribution of 

(x*(0), x*(l), . . . , x*(k), x(k), x{k + l),...,x{k + /)) (2.8) 

coincides with the stochastic evolution of m* under transition matrices 

( P*,...,P* ,A,P,. P). (2.9) 

k I 

Exactly the same statement holds for the Markov chain P& and the ini- 
tial condition m*(x*)A(x*,x) with A replaced by A in the above sequence of 
matrices. 

Proof of Proposition \2.2[ Successive application of the first relations of 
Lemma 12.11 to evaluate the sums over x*(k + I), . . . ,x*(k + 1), and of the 
second relations to evaluate the sums over x(l), . . . , x{k — 1). □ 

Note that Proposition 12.21 also implies that the joint distribution of x*(k) 
and x(k) has the form m* k {x*{k))A{x*{k) 1 x(k)), where ml is the result of 
/c-fold application of P* to m*. 

The above constructions can be generalized to the nonautonomous situ- 
ation. 

Assume that we have a time variable t G Z, and our state spaces as well 
as transition matrices depend on t, which we will indicate as follows: 

S(t), S*(t), P(x,y\t), P*(x*,y*\t), A(x*,x\t), Pit), P*(t), A(t). 

(2.10) 

The commutation relation (1.3) is replaced by A(t)P(t) = P*(t)A(t + 1) or 

A( x \ y \t) : = A ( x *> x I £ ) p ( x > y I = p *( x *>y*\ f ) A ^*> y l 

xes(t) y *es*(t+i) 

(2.11) 
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Further, we set 



s A (t) 
s A (t) 



{(x*,x) G S*(t) x S(t) \A(x*,x\t) > 0}, 
{(x*, x) G S*(t) x S(t + 1) | A(x*, x\t)> 0}, 



(2.12) 



and 



P A ((x*,x),(y*,y)\t) 



P{x,y\t)P*{x*,y* \ t)A(y*,y\t + 1) 
A(x*,y 1 t) 

0, 




A{x*,y\t) >0, 



P A ((^,z),GA?/)|t) 



P(x, y 1 1 + l)P*(x*, | t)A(y*, x\t+l) 
A(x*,x 1 1) 



The nonautonomous generalization of Proposition 12.21 is proved in exactly 
the same way as Proposition 12.21 Let us state it. 

Proposition 2.3. Fix to G Z ; and let m*(x*) be a probability measure on 
S*(t ). Consider the evolution of the measure m(x*)A(x*,x\t ) on S\(to) 
under the Markov chain P\(t), and denote by (x*(t +j),x(t Q +j)) G S^to+j) 
the result after j = 0, 1, 2, . . . steps. Then for any k,l = 0, 1, . . . the joint 
distribution of 



{x*(t ),x*(t +l), x^tQ+^.xito+^.xitQ+k+l), x(t +k+l)) (2.15) 



coincides with the stochastic evolution of m* under transition matrices 



P*(t ), . ..,P*(t + k — 1), A(t + k), P(t + k),..., P(t + k + I - 1) (2.16) 



(for k = I = only A(to) remains in this string). 

A similar statement holds for the Markov chain P&(t) and the initial 
condition m*(x*)A(x* , x \ to): For any k, I — 0, 1, . . . the joint distribution of 



(x*(t ),x*(t + l), . . . ,x*(t + k),x( y t + k + l),x( y t + k + 2), . . . , x(t + k + l + l)) 



coincides with the stochastic evolution of m* under transition matrices 



P*(t ),...,P*(t + k-l),A(t + k),P(t + k + l),...,P(t + k + l). (2.18) 



Remark 2.4. Observe that there is a difference in the sequences of times 
used in (I2.8P and (12. 17ft . The reason is that for nonautonomous Pa, the state 
space at time t is a subset of S*(t) x S(t + 1), and we denote its elements as 
(x*(t),x(t + 1)). In the autonomous case, an element of the state space <Sa 
at time t was denoted as (x*(t), x(t)). 



(2.17) 
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2.2 Multivariate Markov chains 

We now aim at generalizing the constructions of Section 12.11 to more than 
two state spaces. 

Let <Si, . . . , S n be discrete sets, Pi, . . . , P n be stochastic matrices defining 
Markov chains on them, and let Af, ...,A™_ X be stochastic links between 
these sets: 

P k : S k x S k -> [0, 1], ) j Ph{x,y) = 1, x e S fc , fc = l,...,n; 

2/SS fc 

A^_! : <S fe x 5 fe _i -> [0,1], Aj^_ 1 (x, 3 /) = 1, x e <S fe , fc = 2,...,n. 

(2.19) 

Assume that these matrices satisfy the commutation relations 

Ati := AtiPfc-i = PfcAti, k = 2, . . . ,n. (2.20) 
The state spaces for our multivariate Markov chains are defined as follows 



n 

5 i n) = {(^ ■ ■ ■ ,x n ) e Si x ■ ■ ■ x S n \ Y[A-k-i(xk,x k -i) i- o|, 

fc=2 
n 

5 i n) = {(^l, •••>&«) e 5i x • • • x S n I JjA^^arfc.ajfe-O ^ o|. 



(2.21) 



£i=2 



The transition probabilities for the Markov chains P j [ n ' ) and Pj^ are defined 
as (we use the notation X n = (xi, . . . , x n ), Y n = (yi, . . . , y n )) 

{p ( \ rr p k( x k,yk)^- k -i(Vk,yk-i) A Afc / Wn 
^i(^i^i) 11 M 7- - i , 11 ^t 1 (x kl y k - 1 )>0, 
k=2 L± k -i\Xk,yk-l) k=2 
0, otherwise, 

(2.22) 

p£\X n , Y n ) = P{x x , yi ) f[ P *^f-^^ . (2 . 2 3) 

t= 2 ^k-A X k^ X k-l) 

One way to think of pj n) and P| n) is as follows. For P^ , starting from 
X n = (xi, . . . ,x n ), we first choose y x according to the transition matrix 

P(xi,yx), then choose y 2 using ; which is the conditional dis- 

tribution of the middle point in the successive application of P 2 and Af 
provided that we start at x 2 and finish at y x , after that we choose y 3 using 
the conditional distribution of the middle point in the successive application 
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of P3 and A| provided that we start at x 3 and finish at y 2 , and so on. Thus, 
one could say that Y n is obtained by the sequential update. 

For Pj± , starting from X n = (xi, . . . ,x n ) we independently choose 

y 1 ,...,y n according to P^x^yx) for y x and A^fo,,^) > tor V*>, 

k = 2, . . . , n. The latter formula is the conditional distribution of the middle 
point in the successive application of P& and A^_ x provided that we start 
at Xk and finish at Xk-i- Thus, it is natural to say that this Markov chains 
corresponds to the parallel update. 

We aim at proving the following generalization of Proposition 12.21 

Proposition 2.5. Let m n (x n ) be a probability measure on S n . Consider the 
evolution of the measure 

■Al(x 2 ,x 1 ) (2.24) 

on Sj£ under the Markov chain p| , and denote by (xi(j), . . . ,x n (j)) the 
result after j = 0, 1, 2, . . . steps. Then for any k\ > k 2 > ■ ■ • > k n > the 
joint distribution of 

(x n (0) , . . . , , X n —\ {k-a) 1 -En— 1 (J^n -Q j • • • ) •^n— 1 \^n— l) j 

x n - 2 {k n ~i), ■ ■ ■ ,x 2 (k 2 ), x 1 (k 2 ), ■ ■ .fX^ki)) 
coincides with the stochastic evolution of m n under transition matrices 

( P w , . „ , Pn , A"_ 1; P n -l, . . . , P n -i , A"_2, • • • , Aj, Pj, . „ , Pj ). (2.25) 

fc n kn— i A; n fci 

Exactly the same statement holds for the Markov chain P^ and the initial 
condition 

m(x n )A^_ 1 (x n , x n --i) ■ ■ ■ A\{x 2 , xi) (2.26) 
with A 's replaced by A 's in the above sequence of matrices. 

The following lemma is useful. 
Lemma 2.6. Consider the matrix A : S n X i?| n — > [0, 1] gwen by 

A(x n , (xt, . .., Xn-t)) := AJJ_ 1 (a; n , x„_i) ■ • • A^, ffi). (2.27) 
T/ien AP[™ 1 ' ) = P n A. //we denote this matrix by A then 

{ Pj"- 1 ^-!, Y n ^)P n (x n , y n )A(y n , Y^) v 
P£\x n ,Y n ) = l A(x„,F n _i) ' A ^> y ^> U > 

10, otherwise. 

(2.28) 



19 



Also, using the same notation, 

A(x n , X n _i 



P {n) (X Y) = ^ l JJ^zh Yn-l)Pn{Xn,yn)Myn,X n -i) (2 29) 



Proof of Lemma \2.b\ Let us check the commutation relation AP^ n l ' = P n A. 
We have 

AP^ 1) (x n ,y n _ 1 )= ^ ^(Xn.^-O-.-A^Xi) 

Xl,...,X„_l 

where the sum is taken over all xi,...,x n -\ such that 
Yik=2 ^k-i( x k^yk~i) > 0. Computing the sum over xi and using the 
relation AfPi = A^ we obtain 



APjf 1) (x n ,y n _ 1 )= A^_ 1 (x n ,a; n _ 1 )---A^(x3 ) x 2 ) 



X 2 ,...,X n _i 

n-1 



xP2(ar 2j ^)A 1 (^,j/i)|| Afc ^ ~ ^ • (2.31) 

fc=3 



Now we need to compute the sum over x 2 . If A^(x 2 ,?/i) = then 
P 2 (x 2 , y 2 ) = because otherwise the relation A^ = P 2 A^ implies that 
Ai(j/2j3/i) — 0, which contradicts the hypothesis that y n _! G iS>| n 1 . Thus, 
we can extend the sum to all x 2 6 «S 2 , and the relation A|P 2 = A| gives 

Ap| n ~ 1) (x n ,y n _ 1 )= Yl K-i(x n ,x n -i)--*4(x4,x 3 ) 

x P 3 (x 3 , y 3 ) A^s, y^M, V.) ff P *^f^f^ . ( 2 .32) 

Continuing like that we end up with 

K-l(yn~i,y n -2) ■ ■ -Al(y2,yi) ^ A"_i(x n ,x n „i)P n _i(x n _i,?/ n _i), (2.33) 

which, by A™_ 1 P n _ 1 = P n A™_ 1 is exactly P n A(x n ,Y n _i). Let us also note 
that 

A(s„,y n _!) = A^_ 1 (o; n ,2/ n _ 1 )A^:^(y n _i,j/ n _2)---A?(y 2 ,yi). (2.34) 
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The needed formulas for and P^ are now verified by straightforward 
substitution. □ 

Proof of Proposition ^. 51 Let us give the argument for P^] for pj± the proof 
is literally the same. By virtue of Lemma I2.6[ we can apply Proposition 12.21 
by taking 

S = S n , S = iS^ , P = P ni P = Pj^ , = k n , I — k\ — k nj 

(2.35) 

and A(x n ,X n _i) as in Lemma [2761 Proposition 12.21 says that the joint distri- 
bution 

(x n (0),x n (l), . . .,x n (k n ),X n -i(k n ),X n -i(k n + 1), ■ ■ .,X n -i(ki)) (2.36) 
is the evolution of m n under 

(P n ,...,P y A, Pt 1 \...,Pt 1) ) - ( 2 -37) 

k n k\—k n 

Induction on n completes the proof. □ 

As in the previous section, Proposition 12.51 can be also proved in the 
nonautonomous situation. Let us give the necessary definitions. 

We now have a time variable t E Z, and our state spaces as well as 
transition matrices depend on t: 

S k (t), P k (x,y\t), k = l,...,n, A^fa, x k ^ \ t), k = 2,...,n. 

(2.38) 

The commutation relations are 

Ati(t) := Ati(t)JWt) = Pk(t)AU(t + 1), k = 2, . . . , n. (2.39) 
The multivariate state spaces are defined as 

n 

4 n) = {(x l ,...,x n )eS 1 (t)x---xS n (t)\Y[A k k _ 1 (x k ,x k _ 1 \t)^o}, 

k=2 

n 

= {(x 1 ,...,x n ) eStit + n-1) x ••• xS n {t)\ Y[&k-i( x k,Xk-i\t + n 

k=2 

Then the transition matrices for P^ 1 ' and pj± are defined as 
P<?\X n , Y n 1 1) = Pl { Xl , Vl 1 1) Pfc(X " ^ 1 1 + 1} (2.40) 
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if Ylk=2 ^k-i( x ki Vk-i 1 1) > and otherwise; and 
p£ ) (X n ,Y n ) = P(x 1 ,y 1 \t + n-l) 



-A- P k (x k , Vk\t + n- fc)A|_i(yfc, %fc-i 1 t + ra - fc + 1) 



fc=2 



Proposition 2.7. F^r to £ ^? one? /et m n (x n ) be a probability measure on 
S n (to). Consider the evolution of the measure 

m n (x n )A"_ 1 (x n ,x n -i 1 1 ) • • • Al(x 2 ,xi 1 1 ) (2.42) 

on S^\to) under P^ (t). Denote by (xi(t + j), . . . ,x n (t + j)) the result 
after j = 0, 1, 2, . . . steps. Then for any k\ > k 2 > ■ ■ ■ > k n > t the joint 
distribution of 

(x n (to), ■ ■ ■ , X n (k n ), X n —i(k n ), X n -\{k n + 1), . . . , X n _i(/c n _i), 

x n _ 2 (k n -i), ■ ■ • , x 2 {k 2 ),x 1 (k2), • . • , x^h)) 
coincides with the stochastic evolution of m n under transition matrices 

P„(t ), P n {k n - 1), A^M, Pn-l{k n ), P n -l{K-l - 1), 

A^(/c n _i), . . . , A\{k 2 ), P x {k 2 ), . . . , P\{ki - 1). 

A similar statement holds for the Markov chain P^\t) and the initial 
condition 

m{x n )/S^ l _ 1 {x n , x n _i 1 1 ) ■ ■ • &l{x 2 , xt 1 1 + n - 2). (2.43) 
For any k\ > k 2 > ■ ■ ■ > k n > t the joint distribution of 

(x n (t ), . . . ,x n (k n) ; x n —\ {k n -\- 1) , x n ^\(k n -\- 2) , . . . , x n —\ {k n —\), 

x n -2{K-i + 1), • • - ,x 2 {k 2 ),xi{k 2 + 1), . . . ,ari(fci)) 

coincides with the stochastic evolution of m n under transition matrices 

P n (t ), P n {k n - 1), K-l^n), Pn-^kn + 1), . . . , Pn-l{K-l ~ 1), 

K-MK-i), • • • , A?(A; 2 ) ) P^ + 1), . . . , Pxfa - 1). 
The proof is very similar to that of Proposition 12.51 



22 



2.3 Toeplitz-like transition probabilities 

The goal of this section is to provide some general recipe on how to construct 
commuting stochastic matrices. 

Proposition 2.8. Let ct\, . . . , a n be nonzero complex numbers, and let F(x) 
be an analytic function in an annulus A centered at the origin that contains 
all aj 1 's. Assume that F^a^ 1 ) ■ ■ ■ F(a~ l ) ^ 0. Then 



a i%=l 



V 1 1 \ n 1 yi< ... <yn eZ 

(2.44) 

where 

and the integral is taken over any positively oriented simple loop in A. 

Proof of Proposition lKR Since the left-hand side is symmetric with respect to 
permutations of y^s and it vanishes when two y/s are equal, we can extend 
the sum to U 1 and divide the result by n\. We obtain 



2^ det [<*i%=i det - Vi)\lj=i = n! det [ 2^ a lf( x 3 - y)\ k i=1 - 

(2.46) 



yi,—,y n €.'& y=—oo 



Further, 



y= — OQ y= — oo 

1 f +00 y ] P 



2m J w ~ f-^ z x i-v +1 2m J ~ z x i~y +1 

\z\=ci<\a k \-l y ~ X]+1 \z\=c 2 >\a k \-l y ~ °° 



27ri / 1 - a k z 2m J 1 - a k z k k 

|z|=c 1 <|a fc |- 1 |2|=c 2 >|a fc |- 1 

□ 

Proposition 2.9. In the notation of Proposition HOI assume that the vari- 
able y n is virtual, y n = virt, and set f(x k — virt) = a^ k for any k = 1, . . . , n. 
Then 

Ffa-^-^Fia- 1 ) ^ ^ "''-'-^ det " Vi) ^ =1 = det ^G^ 1 " 

V 1 ^ V n_1 ^ !/l<-<tfn-ieZ 

(2.47) 
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Proof of Proposition [Ql Expansion of det [f(xj — yi)]"- =1 along the last row 
gives 

n 

det [f( Xj - = J2(-l) n ~ k < k ■ det [f( Xi - yi )]i=i,..,n-i 

j=l,...,k—l,k+l,...,n 

(2.48) 

The application of Proposition 12.81 to each of the resulting summands in the 
left-hand side of the desired equality produces the expansion of det [cK^] i - =1 
along the last row. □ 

For n = 1, 2, . . . , denote 



X n = {{xi, . . . , x n ) G Z n | Xl < ■ ■ ■ < x n }. (2.49) 

In what follows we assume that the (nonzero) complex parameters 
ai,a 2 ,... are such that the ratios det[a* J ]". =1 / det [a^ -1 ^ . =1 are nonzero 
for all n = 1, 2, . . . and all (xi . . . , x n ) in X n . This holds, for example, when 
all Oj's are positive. The Vandermonde determinant in the denominator is 
needed to make sense of detfo!^^ . =1 when some of the ay's are equal. 

Under this assumption, define the matrices X n x X n and X n x X n -\ by 

T n (a h ...,a n ;F){X,Y) = J — — ' , X, Y G X n , 

det[o;/] iJ=1 Ilj=iF{<Xj ) 

T n _ l (a 1 ,...,a n ,F)(X,Y) - > X G X n , F G X n _ 

where in the second formula y n = virt. By Propositions 12.81 and 12.91 the sums 
of entries of these matrices along rows are equal to 1. We will often omit the 
parameters atj from the notation so that the above matrices will be denoted 
asT„(F) and T^iCF). 

We are interested in these matrices because they have nice commutation 
relations, as the following proposition shows. 

Proposition 2.10. Let Fx and F 2 be two functions holomorphic in an an- 
nulus containing aj 1 's, that are also nonzero at these points. Then 

T n {F x )T n {F 2 ) = T n (F 2 )T n (Fi) = T^F,), 

UFJT^F,) = T^F^T^F,) = T™_ x {FxF 2 ). 1 ' } 

Proof of Proposition WUK The first line and the relation T™_ 1 (F 1 )T n _ 1 (F 2 ) = 
Tn-i{F\F 2 ) are proved by straightforward computations using the fact the 
Fourier transform of F\F 2 is the convolution of those of F\ and F 2 . The 



24 



only additional ingredient in the proof of the relation T n (F 1 )T"_ 1 (F 2 ) = 
^(fiFa) is 

2>(* - v)Mv - virt ) = E ^ - v)< = Fii^x ( 2 - 51 ) 

3/SZ y6Z 

□ 

Remark 2.11. In the same way one proves the commutation relation 

Z^iM^W = ^(F^T^F) (2.52) 
but we will not need it later. 

2.4 Minors of some simple Toeplitz matrices 

The goal of the section is to derive explicit formulas for T n (F) and T^_ X {F) 
from the previous section for some simple functions F. 

Lemma 2.12. Consider F(z) = 1 +pz, that is 

(p, m = 1, 

/(m) = h, m = 0, (2.53) 




T/ien /or integers x\ < • ■ ■ < x n and y 1 < ■ ■ ■ < y n 

Jpl*=i(*i-w) j/^ - Xi 6 {-1,0} /or aW 1 < i < n, 
det [f( Xi - y 3 )} Lj=l = < 

J IU, otherwise. 

(2.54) 

Proof of Lemma Y2.1£\ If Xj < yj for some z then %k < yi for < i and / > i, 
which implies that f(x^ — yi) =0 for such k, I, and thus the determinant in 
question vanishes. If x^ > yi + 1 then x^ > yi + 1 for > z and / < i, which 
means f{xk — yi) = 0, and the determinant vanishes again. Hence, it remains 
to consider the case when X{ — yi G {0, 1} for all 1 < i < n. 

Split {xj}™ =1 into blocks of neighboring integers with distance between 
blocks being at least 2. Then it is easy to see that det \f{xi — yj)] splits into 
the product of determinants corresponding to blocks. Let (xk, ■ ■ ■ , be 
such a block. Then there exists m, k < m < I, such that Xi — y^ + 1 for 
I < i < m, and Xj = |/j for m < i < /. The determinant corresponding to 
this block is the product of determinants of two triangular matrices, one has 
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size m — k and diagonal entries equal to p, while the other one has size I — m 
and diagonal entries equal to 1. Thus, the determinant corresponding to this 
block is equal to p m ~ h , and collecting these factors over all blocks yields the 
result. □ 



Lemma 2.13. Consider F(z) = (1 — qz) 1 , that is 

f(m) 



q, 

0, otherwise. 



(2.55) 



(i) For integers x\ < ■ ■ ■ < x n and y± < ■ • • < y n 



gE£=l(*i w) x < y < x i < i < n 

deb\f[x i -y j )] itj=1 = <j Q , , (2.56) 



otherwise. 



(The condition xq < y\ above is empty.) 

(ii) For integers x\ < ■ ■ ■ < x n and y± < 
variable y n = virt such that f(x — virt) = q x , 



< y n -i, and with virtual 



det [f(xi -yj)] 



l)«-i g ELi^-E" 



: , Xi < yt < x i+1 , 1 < i < n - 1, 
otherwise. 

(2.57) 

Proof of Lemma {MUM (i) Let us first show that the needed inequalities are 
satisfied. Indeed, if Xi < yi for some % then det \f{xi — yj)} = by the same 
reasoning as in the previous lemma. On the other hand, if > yi then 
Xk>yilork>i — l,l<i. Let i be the smallest number such that > y^. 
Then columns % and z + 1 have the form 



qXi-\~yi-i 

qXi-l-Vi 



q Xi yi- 1 * * 



(2.5J 



where the 2x2 block with powers of q is on the main diagonal. This again 
implies that the determinant vanishes. On the other hand, if the interlac- 
ing inequalities are satisfied then the matrix [f(xi — yj)] is triangular, and 
computing the product of its diagonal entries yields the result. 

(ii) The statement follows from (i). Indeed, we just need to multiply 
both sides of (i) by q Vl , denote yi(< x±) by virt, and then cyclically permute 
y/s. □ 
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Lemma 2.14. Consider F(z) = p + qz(l — qz) 1 , that is 



m = 0, 



f(m) 



q m , m > 1, 



(2.59) 



0, otherwise. 

(i) For integral x\ < ■ ■ ■ < x n and y± < ■ ■ ■ < y n 

det [f( Xi - Vj)}l j=1 = gi:?=i(^-w)p#{*[»*=I«}(l _ p )#{i|^-i=2/i> (2.60) 

if Xi-i <Ui<Xi for all 1 < i < n, and otherwise. 

(ii) For integral x\ < ■ • • < x n and y\ < ■ ■ ■ < y n -i, and with virtual 
variable y n = virt such that f(x — virt) = q x , 

det [f(xi - %-)]" i=1 = (-l) n - 1 g Sr= 1 ^-E?= 1 1 i/i p #{i|^ + i=w}( 1 _p)#{*l«*=w> 

(2-61) 

^/^i < 2/i < /or a/Z 1 < i < n — 1, and otherwise. 

Proof of Lemma 2.14 (i) The interlacing conditions are verified by the same 
argument as in the proof of Lemma I2.13( i) (although the conditions them- 
selves are slightly different). Assuming that they are satisfied, we observe 
that the matrix elements of [f(xi — yj)] are zero for j > i + 2 because 
%i < Vi+i < Vi+2 an d fijn) = for m < 0. Further, the + l)-element is 
equal to p if Xi = y^\ or if < y^\. Thus, the matrix is block-diagonal, 
with blocks being either of size 1 with entry f(xi — yi), or of larger size having 
the form 



x k~Vk 



V 





q x k+i—yk g x k+i—yk+i p 
q x k+2—yk q x k+2—yk+i q x k+2—yk+2 



\ q 



X l-Vk 



q 



x l~Vk+l 



q 



x l~Vk+2 









■, x i-yi 



\ 



J 



(2.62) 



with x k = y k+1 , ...,xi-! = yi, and x fc _i < y k , x { < y l+v The determinant of 
( I2.62p is computable via Lemma 1.2 of [6], and it is equal to 



-, x i-yk i 



1-p) 



l-k 



jX k -\ hxi-(y k +-'+yi)n ~ p) 



l-k 



Collecting all the factors yields the desired formula. 

The proof of (ii) is very similar to that of Lemma l2.13( ii). 



(2.63) 



□ 
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2.5 Examples of bivariate Markov chains 

We now use the formulas from the previous two sections to make the con- 
structions of the first two sections more explicit. 

Let us start with bivariate Markov chains. Set S* = X n and S = £ n -i, 
where the sets X m , m = 1,2, ... , were introduced in Section 12.31 We will 
also take 

A = ^(ai, . . . , a n ; (1 - a n z)- 1 ) (2.64) 

for some fixed ai, . . . , a n > 0. 
The first case we consider is 

P = T n _i(ai, . . . ,a n _i; 1 + f3z), P* = T n (ai, . . . , a n ; 1 + f3z), ft > 0. 

(2.65) 

Then Proposition 12.101 implies that 

A = AP = P*A = T"_ 1 (ai, ...,a»; (1 + - a n z)). (2.66) 

According to (12.221) . (I2.23p . we have to compute expressions of the form 

P*(x*,y*)A(y*,y) P*(x*, y*)A(y*, x) 

A(x*,y) A(x*,x) 

for the sequential and parallel updates, respectively. 

We start with the condition probability needed for the Markov chain Pa. 

Proposition 2.15. Assume that x* G S* and y G S are such that 
A(x*,y) > 0, that is, x* k < y^ < x* k+1 for all 1 < k < n — 1. Then the 
probability distribution 

P*(x*,y*)A(y*,y) 

a / ; \ , V e S , 2.67 

A(x*,y) 

has nonzero weights iff 

vl ~ x l e {-1,0}, yk-i<y* k <yk, k = l,...,n, (2.68) 

(equivalently, m&x(x* k — l,yk-i) < y| < mm(x* k ,yk — 1) for all k), and these 
weights are equal to 

max(x fc -l,y fc _ 1 )<mm(x fc ,y fc -l) 
k=l,...,n 

with empty product equal to 1. 
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Remark 2.16. One way to think about the distribution of y* G S* is as fol- 
lows. For each k there are two possibilities for y k : Either max(x£ — 1, y^-i) = 
mm(x* k ,yk — 1), in which case y k is forced to be equal to this number, or 
max(x£ — l,yk-i) = x* k — 1 and mm(x* k ,yk — 1) = x* k , in which case y k is 
allowed to take one of the two values x* k oi x* k — 1. Then in the latter case, 
x* k — y k are i. i. d. Bernoulli random variables with the probability of the 
value equal to a n /(a n + (3). 

Proof of Proposition \2.15\ The conditions for non- vanishing of the weights 
follow from those of Lemmas 12.121 and I2.13[ namely from (12.541) and (12.571) . 
Using these formulas we extract the factors of P*(x* ,y*)A(y* ,y) that depend 
on y*. This yields (ct„//3)^ i;=li/i *. Normalizing these weights so that they 
provide a probability distribution leads to the desired formula. □ 

Let us now look at the conditional distribution involved in the definition 
of the Markov chain The following statement is a direct consequence of 
Proposition 12.151 

Corollary 2.17. Assume that x* G S* andx G S are such that A(x*,x) > 0, 
that is, x* k < Xk < x* k+1 for alll < k < n—1. Then the probability distribution 

P*(x*,y*)A(y*,x) 



A(x*, x) 



y* G S\ (2.70) 



has nonzero weights iff max(x* k — l,Xk-i) < yl < mm(x* k ,Xk — 1), and these 
weights are equal to 

-r-r / 8 \ x * k ~ y * k ( a \ 1 ^ +y ^ 

, • . , . A^h) ■ (2 - n: 

max(Xj- l,x k _ 1 )<mm(xl,x k — 1) 
k=l,...,n 



Let us now proceed to the case 

P = T n _i(«i, . . . , a n _i; (1 - 7^) _1 ), P* = T n (ai, . . . , a n ; (1 - 7^) -1 )- 

(2.72) 

We assume that < 7 < min{a;i, . . . , a n }. 
By Proposition 12.101 

A = AP = P*A = T^at, ...,a n ; 1/((1 - a n z){l - jz))). (2.73) 

Again, let us start with Pa- 

Proposition 2.18. Assume that x* G S* and y G S are such that 
A(x*,y) > 0, that is, xl_ t < y^ — 1 < x* k+1 for all k. Then the probabil- 
ity distribution 

^pMfA^ ,* e5 *, (2 . 74) 

A{x*,y) K J 
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has nonzero weights iff 



4-i < vl < 4> Vk-i<yl<yk, k = i,...,n-i, (2.75) 

(equivalently, max(4_i + l> 2/fc-i) < y| < min(x^, — 1) for all k), and these 
weights are equal to 

f[ r^4^ • (2-76) 

J. J. mm(a'*,y fc -l) 

£ (On/7)' 
«=max(x*_ 1 +l,j/ fc _ 1 ) 

i/ere max(iQ + is assumed to denote — oo. 

Remark 2.19. Less formally, these formulas state the following: Each y\ 
has to belong to the segment [max(i^_ 1 + 1, yk-i), min(4> Vk — 1)], and the 
restriction that A(x*,y) > guarantees that these segments are nonempty. 
Then the claim is that yl's are independent, and the distribution of y* k in 
the corresponding segment is proportional to the weights (a n / r y) Vk . In other 
words, this is the geometric distribution with ratio a„/7 conditioned to live 
in the prescribed segment. 

Proof of Proposition \2.18l Similarly to the proof of Proposition 12.151 we use 
Lemmas 12.131 to derive the needed inequalities and to single out the part of 
the ratio P*(x*, y*)A(y*, y)/A(x*, y) that depends on y*. One readily sees 
that it is equal to {ot n /l) =1 , and this concludes the proof. □ 

Let us state what this computation means in terms of the conditional 
distribution used in the construction of Pa- 

Corollary 2.20. Assume that x* £ S* andx £ S are such that A(x*,x) > 0, 
that is, 4-i < Xk — 1 < 4+1 f or a ^ k- Then the probability distribution 

P, (*;f'> A f*>, y - e5 -, (2.77) 
A(x*,x) y V ' 

has nonzero weights iff ma.x(x k \_ 1 + 1, Xk-i) < y* k < min(4> x^ — 1) for all k, 
and these weights are equal to 

TT r^P • (2-78) 

J. J. m'm^^.Xh — 1) 

£ K/t)' 

i=max(a:^_ 1 +l,xj._ 1 ) 
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In the four statements above we computed the ingredients needed for 
the constructions of the bivariate Markov chains for the simplest possible 
Toeplitz-like transition matrices. In these examples we always had x* k > yl, 
or, informally speaking, "particles jump to the left". Because of the previ- 
ous works on the subject, it is more convenient to deal with the case when 
particles "jump to the right". The arguments are very similar, so let us just 
state the results. 

Consider 

P = T n - 1 (a 1 ,...,a n - 1 ;l + /3z- 1 ), P* = T n (a u . . . , a n ; 1 + /J*" 1 ), (3 > 0. 

(2.79) 

• For Pa, we have max(^, Dk-i) < ul < min(a;£ + 1,^ — 1). This segment 
consists of either 1 or 2 points, in the latter case y* k — x* k are i. i. d. Bernoulli 
random variables with the probability of equal to (1 + an/3) -1 . 

• For Pa, we have max(x£, Xk-i) < yl < mm(x* k + 1, Xk — 1), and the rest 
is the same as for Pa. 

Now consider 

P = T n _i(ai, . . . , a n _i; (1 - 7^T 1 )~ 1 ), P* = T n (a u . . . , a n ; (1 - 7^r 1 )~ 1 ), 

(2.80) 

for < 7 < minjo;^ 1 , . . . , a^ 1 }. 

• For Pa, we have max^,^!) < yl < min(x* k+1 , yk) — 1, and yl are in- 
dependent geometrically distributed with ratio (a n j) random variables con- 
ditioned to stay in these segments. 

• For Pa, we have m&x(x* k , Xk-i) < yl < min(x£ +1 , Xk) — 1, and the rest is 
the same as for Pa. 

Thus, we have so far considered eight bivariate Markov chains. It is 
natural to denote them as 

P A (l+/3z ±1 ), P a (1+^ ±1 ), P a ((1-7^ 1 )' 1 ), Pa((1-7^ ±1 )- 1 ). 

(2.81) 

Observe that although all four chains of type Pa live on one and the same 
state space, all four chains of type Pa live on different state spaces. For the 
sake of completeness, let us list those state spaces: 

5a = { (x* , x) G X n x I x* k + 1 < Xk < x* k+1 for all k} 

Sa(1 + (3z) = {(x*, x) G X n x I x* k < Xk < x*k+i f° r a ^ ^} 

Sa(1 + f3z~ l ) = { (x* , x) G X n x 3£ n _i I x* k + 1 < Xk < x* k+1 + 1 for all k} 

Sa((1 — l z )~ l ) — { (x* , a;) G X n x 3£ n _i | x* k _ 1 + 2 < Xk < x* k+1 for all k} 

5a((1 — 7-2 -1 ) -1 ) = { (x* , x) G X n x I x£ + 1 < x k < x* k+2 — 1 for all k} 

In the above formulas we always use the convention that if an inequality 
involves a nonexistent variable (like xq or it is omitted. 
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2.6 Examples of multivariate Markov chains 

Let us now use some of the examples of the bivariate Markov chains from 
the previous section to construct explicit examples of multivariate (not nec- 
essarily autonomous) Markov chains following the recipe of Section 12.11 

For any m > we set S m = 3i m , which is the set of strictly increas- 
ing m-tuples of integers. In this section we will denote these integers by 
x™ < ■ ■ ■ < x™. 

Fix an integer n > 1, and choose n positive real numbers a\, . . . , a n . We 
take the maps A^_ x to be 

Ati = TtMi, ■ ■ ■ , a k ] (1 - ockz)- 1 ), k = 2, . . . , n. (2.82) 

We consider the Markov chain i.e., the sequential update, first. Its 
state space has the form 

n 

= {(x\...,x n )eS l x---xS n \ H A^_ x (x m , x m ~ l ) > o} 

m=2 (2.83) 

= \{x%} m =i,...,n C | x™ +1 < x% < x™+i for all k,m\. 

< fe=l,...,m > 

In other words, this is the space of n interlacing integer sequences of length 
l,...,n. 

Let t be an integer time variable. We now need to choose the transition 
probabilities P m (t), m = 1, . . . , n. 

Let {F t (z)} t >t be a sequence of functions each of which has one of the 
four possibilities: 

F t (z) = (l + fl+z) or (l + Pt/z) or (1- 7 ^) _1 or (1-^/z)- 1 . (2.84) 
Here we assume that 

Pf, it > °> it < min{ai, . . . , a n }, % < min^ 1 , . . . , a" 1 }. (2.85) 
We set 

P m (t) = T m (a 1 ,...,a m ;F t (z)), m = l,...,n. (2.86) 

Then all needed commutation relations are satisfied, thanks to Proposi- 
tion EM 

The results of Section [2~5l enable us to describe the resulting Markov chain 
on <S^ as follows. 

At time moment t we observe a (random) point {x™(t)} G . In order 
to obtain {x™(t + 1)}, we perform the sequential update from level 1 to level 
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n. When we are at level m, 1 < m < n, the new positions of the particles 
x™ < ■ ■ ■ < x™ are decided independently. 

(1) For F t (z) = 1 + /3fz, the particle x™ is either forced to stay where 
it is if x™Z\(t + 1) = x™(t), or it is forced to jump to the left by 1 if 
X T~ (^ Jr 1) = x T(t), or it chooses between staying put or jumping to the 
left by 1 with probability of staying equal to 1/(1 + /^a" 1 ). This follows 
from Proposition 12.151 

(2) For F t (z) = 1 + fit I 'z, the particle x™ is either forced to stay where 
it is if x™~ x (t + 1) = x™(t) + 1, or it is forced to jump to the right by 1 if 
x ^-i {t+ 1) = x™(t) + 1) or ft chooses between staying put or jumping to the 
right by 1 with probability of staying equal to 1/(1 + (3fa m ). 

(3) For F t (z) = (1 — •y^z)~ 1 , the particle x™ chooses its new position 
according to a geometric random variable with ratio a m /jf conditioned to 
stay in the segment 

[max^W + l,xri(t + l)),min(xT(t),x^ l (t + 1) - 1)]. (2.87) 

In other words, it tries to jump to the left using the geometric distribution 
of jump length, but it is conditioned not to overcome x ^Li(t) + 1 (in order 
not to "interact" with the jump of x™_^), and it is also conditioned to obey 
the interlacing inequalities with the updated particles on level m — 1. This 
follows from Proposition 12.181 

(4) For F t (z) = (1 — 7^/^) _1 , the particle x™ chooses its new position 
according to a geometric random variable with ratio a. m lt conditioned to 
stay in the segment 

[max(xr(t),^- 1 (*+l)),min(xr +1 (t),xr 1 (t + l))-l]. (2.88) 

In other words, it tries to jump to the right using the geometric distribution 
of jump length, but it is conditioned not to overcome x™ +1 (t) — 1 (so that it 
does not interact with jumps of x™ +1 ), and it is also conditioned to obey the 
interlacing inequalities with the updated particles on level m — 1. 

Projection to {a;™} m >i. A remarkable property of the Markov chain 
P| n ' ) with steps of the first three types is that its projection onto the 
n-dimensional subspace {x\ > x\ > ■ ■ ■ > x™} (the smallest coordinates 
on each level) is also a Markov chain. Moreover, since these are the leftmost 
particles on each level, they have no interlacing condition on their left to be 
satisfied, which makes the evolution simpler. Let us describe these Markov 
chains. 

At time moment t we observe {x\(t) > x\(t) > ■ ■ ■ > x™(t)}. In order to 
obtain {x™(t + 1)}^ =1 , we perform the sequential update from x\ to x™. 
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(1) For F t (z) = 1 + (3fz, the particle x™ is either forced to jump (it is 
being pushed) to the left by 1 if x™ _1 (t + 1) = x™(t), or it chooses between 
not moving at all or jumping to the left by 1 with probability of not moving 
equal to 1/(1 + /^a" 1 ). 

(2) For Ft(z) = 1 + (3f /z, the particle x™ is either forced to stay where 
it is if x™ _1 (t + 1) = x™(t) + 1, or it chooses between staying put or jumping 
to the right by 1 with probability of staying equal to 1/(1 + /3fa m )- 

(3) For F t (z) = (1 — jfz)' 1 , the particle x™ chooses its new position 
according to a geometrically distributed with ratio 7 t + /« n jump to the left 
from the point min(acf (t), x^(t + 1) - 1). That is, if xf{t) < x™~\t + 1) 
then x™ simply jumps to the left with the geometric distribution of the jump, 
while if x?(t) > x^ x {t + l) then x™ is first being pushed to the position 
x™~ (t + 1) — 1 and then it jumps to the left using the geometric distribution. 

(4) For the transition probability with F t (z) = (l—j^/z)" 1 , the particle 
x™ is conditioned to stay below min(x™(£), x™ _1 (t + 1)) — 1, which involves 
x™, thus the projection is not Markovian. 

The Markov chains on {x\ > ■ ■ ■ > x n } corresponding to 1 + (3fz and 
1 + flt / z are the "Bernoulli jumps with pushing" and "Bernoulli jumps with 
blocking" chains discussed in [20]. 

Projection to {x™} m >i. Similarly, the projection of the "big" Markov 
chain to {x\ < x\ < ■ ■ ■ < x™} is Markovian for the steps of types 
one, two, and four, but it is not Markovian for the step of the third type 
F t (z) = (l- 7 +z)- 1 . 

Let us now consider the parallel update Markov chain , or rather one 
of them. 

Choose a sequence of functions G t (z) = 1 + (3 t z~ l with (3 t > 0, and set 

P m (t) = T m (at, . . . ,a m ;G t (z)), m = l,...,n. (2.89) 

In case (3 t = 0, P m (i) is the identity matrix. As before, the needed commu- 
tation relations are satisfied by Proposition 12.101 

The (time-dependent) state space of our Markov chain is 

n 

St\t) = {(x 1 ,...,x n )eS l x---xS n \ Yl AI-i(^, x m ~ l I t+n-m) > o} 

m=2 

f m n(n+l) m m—l m 

= \ { x k }m=l,...,n C Z 2 | X k < X k < x k+l ^ Pt+n-m = 0, 
^ k=l,...,m 

^<^ 1 <^+i + l if Pt+n-m>0}. (2.90) 
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The update rule follows from the analog of Corollary 12 . 1 71 for (1 + (3 t z r ). 
Namely, assume we have {x™(t)} G S^\t). Then we choose {x™(t)} inde- 
pendently of each other as follows. We have 



mBx(xf(t),x^(t)) < xf(t + 1) < inin(sj»(f) + l,x™~\t) - 1). (2.91) 



This segment consists of either 1 or 2 points, and in the latter case 
x™ +1 (t + 1) has probability of not moving equal to (1 + a m /?f_ n+m ) _1 , and 
it jumps to the right by 1 with remaining probability. In particular, if 
(3t-n+m = then x™(t + 1) = x%(t) for all k = 1, . . . , m. 

Less formally, each particle x™ either stays put or moves to the right 
by 1. It is forced to stay put if x™(t) = x™ _1 (t) — 1, and it is forced to 
move by 1 if x™(t) = x™Z\{t) — 1- Otherwise, it jumps with probability 

1 - (1 + anA-n+m)"" 1 - 

Projection to {x™} m >i- Once again, the projection of this Markov chain 
to {x\ > ■ ■ ■ > x™} is also a Markov chain, and its transition probabilities 
follows: Each particle x™ at time moment t is either forced to stay if 
x™{t) = x™~ 1 (t) — 1 or it stays with probability (1 + ctnPt-n+m)^ 1 and jumps 
to the right by 1 with complementary probability. This Markov chain has no 
pushing because x^'s do not have neighbors on the left. This is the "TASEP 
with parallel update", see e.g. [10]. 

Projection to {a^} m >i. We can also restrict our "big" Markov chain 
to the particles {x\, x%, ■ ■ ■ Then at time moment t they satisfy the 

inequalities 



xZz{(t) < x2(t) if Pt+n- m = 0, ClW<<W+l if A+n-m>0, 



and the update rule is as follows. If x™_{(t) = x™(t) + 1 then a;™ moves to 
the right by 1: x™(t + 1) = x™(t). However, if x%z{(t) < x™(t) then 
stays put with probability (1 + a n /5 4 _ n+m ) _1 , and it jumps to the right by 1 
with the complementary probability. 
In the special case when all a, = 1, 



and with the densely packed initial condition x™(n — m) = k — m — 1, the 
Markov chain discussed above is equivalent to the so-called shuffling 
algorithm on domino tilings of the Aztec diamonds that at time n produces 
a random domino tiling of the diamond of size n distributed according to the 
measure that assigns to a tiling the weight proportional to (3 raised to the 
number of vertical tiles, see [39] . 



(2.92) 




(2.93) 
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2.7 Continuous time multivariate Markov chain 



The (discrete time) Markov chains considered above admit degenerations 
to continuous time Markov chains. Let us work out one of the simplest 
examples. 

As in the previous sections, we fix an integer n > 1 and n positive real 
numbers aci, . . . , a n , and take 

Ati = TliK . . . , a k - (1 - a k z)-% k = 2,...,n. (2.94) 
We will consider a limit of the Markov chain , so our state space is 
SjW = \{x^} m =i,...,n C | x™ +1 <x%< a^+i 1 for all k, m\. (2.95) 

l. fc=l,...,m ' J 

In the notation of the previous section, let us take F t (z) = 1 + (3~ / z for 
a fixed /3- > and i = 1,2,.... Thus, we obtain an autonomous Markov 
chain on whose transition probabilities are determined by the following 
recipe. 

In order to obtain {x™(t + 1)} from we perform the sequential 

update from level 1 to level n. When we are at level m, 1 < m < n, for each 
k = 1, . . . ,m the particle x™ is either forced to stay if x™ _1 (t + l) = 2™(t) + l, 
or it is forced to jump to the right by 1 if x™Z\ {t + 1) = x™(t) + 1, or it 
chooses between staying put or jumping to the right by 1 with probability 
of staying equal to (1 + /5~a m ) _1 . Note that, since particles can only move 
to the right, it is easy to order the elements of the state space so that the 
matrix of transition probabilities is triangular. 

We are now interested in taking the limit j3~ — > 0. 

Lemma 2.21. Let A(e) be a (possibly infinite) triangular matrix, whose ma- 
trix elements are polynomials in an indeterminate e > 0: 

A(e) = A + eAx + e 2 A 2 + . . . , (2.96) 

and assume that Aq = 1. Then for any reR, 

\im(A(e)) [r/€] = exp(rAi). (2.97) 

Proof of Lemma \2.21\ For the finite size matrix the claim is standard, and 
the triangularity assumption reduces the computation of any fixed matrix 
element of (A(e)Y T ^^ to the finite matrix case. □ 

This lemma immediately implies that the transition probabilities of the 
Markov chain described above converge, in the limit j3~ — >• and time rescal- 
ing by f3~, to those of the continuous time Markov chain on whose gen- 
erator is the linear in f3~ term of the generator of the discrete time Markov 
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chain. Let denote this linear term by . Its off-diagonal entries are not 
hard to compute: 



( K}m=l,,n, = 1 (2-98) 

\ fe=l,...,m fc=l,...,m / 



if there exists 1 < a < b, 1 < b < n, < c < n — b such that 



rJ> — r b+1 — ... — T fc+C — T 



and xjj 1 = y™ for all other values of (k, m), and 



L {n) ({xT} m =l,..,n, 
\ k=l,...,m 



771=1, 
fc=l,. 



(2.99) 



in all other cases. 

Less formally, this continuous time Markov chain can be described as 
follows. Each of the particles x™ has its own exponential clock, all clocks 
are independent. When x^-clock rings, the particle checks if its jump by one 
to the right would violate the interlacing condition. If no violation happens, 
that is, if 

x h a <x h - x -\ and x» a <x»+\, (2.100) 

then this jump takes place. If x b a = x h ~ l — 1 then the jump is blocked. On 
the other hand, if x h a = then we find the longest string x h a 



r 6+l 



• • • = x a +c an d move all the particles in this string to the right by one. One 
could think that the particle x b a has pushed the whole string. 
We denote this continuous time Markov chain by V^ n ' . 

(n) 

Similarly to P^ , each of the Markov chains P m on S m also has a con- 
tinuous limit as (3~ — > 0. Indeed, the transition probabilities of the Markov 
chain generated by T m (a±, . . . , a m ; 1 + (3~ / z) converge to (x m , y m G S m ) 

V/p-i i {x n ,y m ) 



(2.101) 



( lim (T m (a 1 ,...,a m ;l + (3 / z)) 

_ det [af }™ =1 det [t*-°?1(v? - xf > 0)/{y? - x?)\]™ =1 



det[«^]„ =1 exp(rE i=1 % 



H,j=l 

Thus, the limit of P m is the Doob /i-transform of m independent Poisson 
processes by the harmonic function h(xi, . . . ,x m ) = det - =1 , cf. [40]. 
Let us denote this continuous time Markov chain by V m , and the above 
matrix of its transition probabilities over time r by V m {r). 

Taking the same limit f3~ in Proposition 12.51 leads to the following state- 
ment. 
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Proposition 2.22. Let m n (x n ) be a probability measure on S n . Consider 
the evolution of the measure 

m n (x n )kl^(x n , x n ~ l ) ■ ■ ■ Al{x 2 , x 1 ) (2.102) 

on under the Markov chain V^ n \ and denote by (x 1 (t) , . . . , x n (t)) the 
result after time t > 0. Then for any 

n — /° < . . . < t c{n) —t° < ■■■ < - t° < ■■■ < t c{2) -t° < ■■■ <t a 

U — l n — l n-l — — l c(n-l) ~ l n-2 — ^ l 2 ~ l l ^ l l 

(2.103) 

(here c(l), . . . , c(n) are arbitrary nonnegative integers) the joint distribution 
of 

x (t°),...,x (tj (t°_ 1 ),a; (t n _ 1 ), . . . , x "*), 

(^n-2)) • • • ) x {^2 (£?)> • • • , a; (t/ ^) 
coincides with the stochastic evolution of m n under transition matrices 

T> (f 1 —f°) V it c(:ri) — / c (™)~ 1 '| \ n 

1 n\ L n l n)i ■ ■ ■ 1 'n\ l n l n ) i 1 *-n-li 

T> (+ 1 -+° ) T> (+ c ( n - 1 ) _ y-cO- 1 )- 1 ^ An-1 

' re-lV^n-l Vlif-i^n-HVl t 'n-l J > Yi n-2' • • • ) 

...,A^ 1 (t}-t?),...,P 1 (tf ) -tf ) - 1 ). 

Remark 2.23. It is not hard to see that if in the construction of pj^ we 
used F t (z) = (1 — j~ / z)^ 1 and took the limit 7" — > then the resulting 
continuous Markov chains would have been exactly the same. On the other 
hand, if we used F t (z) = (1 + (3 + z) or F t (z) = (1 — / -f + z)~ 1 then the limiting 
continuous Markov chain would have been similar to , but with particles 
jumping to the left. 

It is slightly technically harder to establish the convergence of Markov 
chains with alternating steps, for example, 

F 2s (z) = 1 + P + (s)z, F 2s+1 = 1 + p-(s)/z, (2.104) 

because the transition matrix is no longer triangular (particles jump in both 
directions). It is possible to prove, however, the following fact: 

For any two continuous functions a(r) and 6(r) on R + with 
a(0) = 6(0) = 0, consider the limit as e — > of the Markov chain P^ with 
alternating F t 's as above, 

P~(s) = ea(es), (3 + (s) = eb(es), (2.105) 
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and the time rescaled by e. Then this Markov chain converges to a continuous 
time Markov chain, whose generator at time r is equal to a(r) times the 
generator of plus b(r) times the generator of the Markov chain similar 
to but with particles jumping to the left. 

The statement of Proposition 12.221 also remains true, but in the definition 
of the Markov chains V m one needs to replace the Poisson process by the 
one-dimensional process whose generator is a(r) times the generator of the 
Poisson process plus b(r) times the generator of the Poisson process jumping 
to the left. 

2.8 Determinantal structure of the correlation func- 
tions 

The goal of this section is to compute certain averages often called correlation 
functions for the Markov chains and P^ with F t (z) = (1 + (3fz ±l ) or 
(1 — 7 t :t 2; ±1 )~ 1 , and their continuous time counterpart P™, starting from a 
certain specific initial condition. 

As usual, we begin with P^ . The initial condition that we will use is 
natural to call densely packed initial condition. It is defined by 

x£(Q)=k-m-l, k = l,...,m, m = l,...,n. (2.106) 

Definition 2.24. For any M > 1, pick M points 

Xj — (yj,mjitj) G Zx {1, . . . ,n}xZ> or Z x {1, . . . , n} x R> , (2.107) 

j = 1,...,M. The value of the Mth correlation function pu of P^ (or 
Pj ) at (xi, . . . , Km) is defined as 

Pm(xi, • • • , %) = Prob{Por each j = 1, . . . , M there exists a kj, 

1 < kj < uij, such that x™ 1 (tj) = yj}. (2.108) 

The goal of this section is to partially evaluate the correlation functions 
corresponding to the densely packed initial condition. 

Introduce a partial order on pairs (m, t) e {1, . . . , n} x Z> or {1, . . . , n} x 
R>o via 

(mi,ti) -< (m 2 ,t 2 ) iff mi < m 2 , h > t 2 and (mi,ti) ^ (m 2 ,t 2 ). (2.109) 

In what follows we use positive numbers a%, . . . , a n that specify the links 
A^_ 1 as in Section [2T6l and as before we assume that 

Pti it > °> It < min{ai, . . . , a n }, % < mm{ai\ . . . , a" 1 }. 

(2.110) 
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Theorem 2.25. Consider the Markov chain P^ with the densely packed 
initial condition and F t {z) = (1 + (3fz ±l ) or (1 — 7 t ± 2 ±1 ) _1 . Assume that 
triplets Kj = (yj,mj,tj) , j = 1, . . . , M, are such that any two distinct pairs 
(rrij,tj), (m,ji,tj>) are comparable with respect to -<. Then 



p A /(xi, • • • , %) = det [JC(xi, Xj)}f )j=v (2-111) 
where 

r( *\ 1 / dw Uttl Ft(w) 

K( yi , mi , tr, y 2 , m 2 , t 2 ) = -— ^ nnni+i(1 _ Qjw) Ww***)] 



i/ie contours Tq, r a -i are closed and positively oriented, and they include the 
poles and {af , . . . , a" 1 }, respectively, and no other poles. 

This statement obviously implies 

Corollary 2.26. For the Markov chain V^ 1 ' , wi/i iae notation of Theo- 
rem \2.25\ and densely packed initial condition, the correlation functions are 
given by the same determinantal formula with the kernel 

If dw e^-t,)/*, 
K{y u m U T l]y2 ,m 2 ,T 2 ) = -— <b — — = _ l[( mi ,t lH (m 2 , ta )] 



(27ri) 2 / ro Jr a _ 1 e*»A n^U 1 - M ^ 2+1 ™ - 2 ' 

Remark 2.27. For the more general continuous time Markov chain described 
in Remark 12.231 a similar to Corollary 12.261 result holds true, where one needs 
to replace the function e t/w by e a(t)/w+b(t)w . 

Proof of Theorem \2.25[ The starting point is Proposition 12.71 The 
densely packed initial condition is a measure on of the form 

m n (x n )A™_ l (x n , x n ~ v ) • • ■ A 2 ^ 2 ,^ 1 ) with m n being the delta-measure at the 
point (— n, —n + 1, . . . , —1) G S n . 

This delta-measure can be rewritten (up to a constant) as 
det[a i i ]ij = i,... >fl det[*^_,(z]J)] fci{= i,... in with 

K-i(x) = ^-J f[(l- aj w)w x+l — , l = l,...,n. (2.112) 

7X1 ^ r o j=l+l " 
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Indeed, Span(\I/"_^ 1 1 = 1, . . . , n) is exactly the space of all functions on Z 
supported by {—1, . . . , — n}. 

We are then in a position to apply Theorem 4.2 of [7]. For convenience of 
the reader, this theorem can be found in Appendix [A] (In fact, the change of 
notation that facilitates the application was already used in Proposition 12.221 
above.) The computation of the matrix M -1 of that theorem follows verba- 
tim the computation in the proof of Theorem 3.2 of [13], where 9j of [13] have 
to be replaced by aj 1 for all j = 1, . . . , n. Arguing exactly as in that proof 
we arrive at the desired integral representation for the correlation kernel. □ 

Finally, one can also derive similar formulas for the Markov chain Pj± . 
As the state space is now 

S%\t) = {{x n {t),x n -\t + 1), . . . ,x\t + n - 1)}, (2.113) 

we need to define the densely packed initial condition differently, cf. the end 
of Section ES We set 

x™(n — m) = k — m — 1, k = 1, . . . , m, m — 1, . . . , n, (2.114) 

and assume that F t {z) = 1 for t = 0, . . . , n — 2. This means that 

A™_ 1 (x m ,x m " 1 |n-m) = A^.^x" 1 , x" 1 " 1 ), m = 2,...,n, (2.115) 

and our initial condition is of the form (12. 43ft . 

Corollary 2.28. For the Markov chain Pj±\ with the above assumptions, 
notation of Theorem \2.25l and densely packed initial condition, under the 
additional assumption that for any two pairs (mj,tj) -< (mjt,tjt) we have 

tj — tji > iriji — rrtj, (2.116) 

the correlation functions are given by the same determinantal formula as in 
Theorem \2.25l 

Proof of Corollary {272^ Comparing the formulas for the joint distributions 
for P^ and P^ in Proposition 12.71 we see that with the densely packed 
initial conditions they simply coincide. Hence, the correlation functions are 
the same. □ 

Note that according to the remark at the end of Section 12.61 the correla- 
tion functions for the shuffling algorithm of domino tilings of Aztec diamonds 
can be obtained from Theorem 12.251 and Corollary 12.281 
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3 Geometry 



3.1 Macroscopic behavior, limit shape 

It is more convenient for us to slightly modify the definition of the height 
function (11.21) by assuming that its first argument varies over Z, and 

h(x,n,t) = \{k\x n k (t)>x}\. (3.1) 

Clearly, this modification has no effect on asymptotic statements. 

We are interested in large time behavior of the interface. The macroscopic 
choice of variables is 

x = [(u-r])L}, n=[r]L], t = tL, (3.2) 

where (z/, rj, r) G K.1 and L ^> 1 is a large parameter setting the macroscopic 
scale. For fixed r] and r, h(x,n,t) = n for v small enough (e.g., v = 0) and 
h(x, n,t) =0 for v large enough. Define the x-density of our system as the 
local average number of particles on unit length in the x-direction. Then, 
for large L, one expects that —L^ 1 dh/du ~ s-density. Thus, our model has 
facets when the x-density is constant (equal to or 1 in our situation), which 
are interpolated by curved pieces of the surface, see Figure [TTBl 

Claim 3.1. The domain T> C where the x-density of our system is 
asymptotically strictly between and 1 is given by 

\Vr- y/rj\ < V"< VV- (3-3) 

Equivalently, x-density 6 (0, 1) iff there exists a (non-degenerate) triangle 
with sides y/v, y/rj, y/r. Denote by n u , ir v and tt t the angles of this triangle 
as indicated in Figure 13.11 Claim 13.11 follows from Proposition 13.21 below. 

The condition (13.31) is also equivalent to saying that the circle centered at 
of radius \frjjr has two disjoint intersections with the circle centered at 1 
of radius \J v/t. In that case, the two intersections are complex conjugate. 
Denote by Q(u,i],t) the intersection in 

e = {^GC|Im(^) >0}. (3.4) 

Then, we have the following properties 

\tt\ 2 = -, \l-n\ 2 = -, arg(fi)=7iv, arg(l - fi) = -7r„. (3.5) 
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o 1 

Figure 3.1: The triangle of ( 13.31) on the left and its scaled version defined by 
intersection of circles on the right. 



The cosine rule gives the angles 7r*'s in (0, n) by 

V + 1] — v 
7i v = arccos | — — 



tt v = arccos I — — — I , (3.o) 



7] + V — 7] 

ir T = arccos | — — 

Proposition 3.2 (Bulk scaling limit). For any k = 1, 2, ... , consider 

= (xj(L),rij(L),tj(L)), j = l,...,k, (3.7) 

such that for any i ^ j and any L > either (rii(L),ti(L)) -< (rij(L),tj(L)) or 
(rij(L),tj(L)) -< (rii(L),ti(L)) (the notation -< was defined in U.3\) ). Assume 
that 

lim — - = v, lim -p- = r), lim -4- = r, j = 1, . . . , k; (3.8) 

L— >oo Li L—*oo L L—>oo Li 

we have (v, T], t) G T>; and also all the differences Xi — Xj, rij — rij, £j — tj do 
not depend on the large parameter L. Then the k-point correlation function 
Pk{xi, • • • , Xfc) converges to the determinant det[i^ ulk ]i<jj<fc ; where 



i^ ik = / d«P * , (3.9) 

where for (rii,ti) -/( i n j^j) the integration contour crosses R + , while for 
(rij, tj) -< (n^, t-,) t/ie contour crosses R_. On t/ie ot/ier hand, if (u, rj, r) ^ P ; 
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then 

lim p k (>ci, ■ ■ ■ , Xk) = 0, if \/v > y/rj + yfr 

lim Pfc(x 1 ,...,x fc ) = 0, if y/v < y/r - y/rj (3.10) 

lim p fc (xi,...,x fc ) = 1, if Vv < y/rj ~ Vt- 

L— >oo 

Proof of Proposition I ff . £l One follows exactly the same steps as in Section 3.2 
of [41], replacing the double integral (35) in there by (14.21) . The deformed 
paths are then like in Figure E31 but with z c = w c . The degenerate cases when 
Xj ^ X> are treated in the same way with limiting kernel K^j being either 
(no residue in the contour integral computation) or triangular (K^j = for 
Xi < Xj) with Ki t i = 1, when the integral in ( 13. 9ft is over a complete circle 
around the origin. □ 



Corollary 3.3. Let p denote the asymptotic x-density. Then, in V, it is 
given by 

P {u,ti,t) = lim Pl {[uL), [r]L),rL) = rrjn e [0, 1]. (3.11) 

L— >oo 

Consequently, 

h( x r m([(u-r,)L],[r,L],TL) 1 /-^+^) 2 

ft(i/, 77, r) := hm = - / tc^v , rj, T)du . 

L^oo L 7T J u 

(3.12) 

Below we perform the integral in (I3.12p to get an explicit expression 
for the limit shape h. Along the way we derive some interesting geometric 
relations. First of all, h is homogeneous of degree one (since it is the scaling 
limit under same scaling in all directions). 

Lemma 3.4. For any a > 0, 

h(au, an, ar) = ah(u, 77, r), (3.13) 

from which it follows 

d d d\ 
v— + r]— + t— h(u, r), t) = h(v, r), r). (3.14) 



dv drj dr 

Proof of Lemma \3.4\ It follows directly from the geometric property 

TTjj^QCIS, CUT], CUT ) = 7T 7 j(l / , T),t). □ 

Therefore, we just need to compute the partial derivatives of h, then the 
limit shape h will be determined by the l.h.s. of ( 13.141) . 
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Proposition 3.5. The partial derivatives of the limit shape h are given by 

dh _ jix^ dh _ ^ _ tx^ dh _ sin(7iy) sin (7^) 
dv ix ' dr] ix ' dr 7rsin(7r T ) 

Remark 3.6. Another expression for the growth velocity is 

dh 1 , 
— = -Imfi {u, v ,r . 3.16 
or 7r 

This can be understood using Proposition 13. 21 The macroscopic growth 
velocity is equal to the average flow of particles, J. It is computed in Sec- 
tion El see with Q = 0: E(J) = -d 2 K(x,n,t;x,n,t). Then, by CTJ]) 
we have E(J) = K(x,n,t;x + l,n,t). Finally, by Proposition 13.21 one gets 
E(J) = ImO/7r. 

As a corollary of Lemma I3~il and Proposition 13.51 the limit shape is given 
as follows. 

Corollary 3.7. For (u,t],t) G T>, we have 

1 / . . sin(7Tj,) sin(7r„) \ 

IX \ sm(7r r ) / 

Proof of Proposition \3.5\ From (I3.12p we immediately have the first relation: 
dh/dv = —tx v /tx. In the derivatives of h with respect to r and rj we have one 
term coming from the boundary term and one from the internal derivative. 
The boundary terms will actually be zero, since the density at the upper 
edge is zero. We need to compute 

d *v = 1 = v-V-t 

dr} ^ArfT- {v - i) -r) 2 ' dr 2t^At]t - {v - rj - r) 2 ' 

Then, we apply the indefinite integrals 

f n dX 5 = arcsin(x/|a|) +C, f = -Va 2 - x 2 + C. (3.19) 

For the derivative with respect to n, 



TX 



fri + r — v\ ( 7] + r — u\ 
= tx 2 + arcsin = ix — arccos ,(3.20) 
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Figure 3.2: (a) Limiting density of the particles with r = 1. (b) The associ- 
ated limiting height function. Two facets are visible. 



the latter being 7r„. Finally, 



7T 



dh 



dr 



l du > + (1 + v ^77) 7r?? ((^ + v ^ ) 2 ) ^ r ) 



a/ 4?]t — [y — r\ — r) 1 
2^ 



r sin 71% 



(3.21; 



and by the sinus theorem for the triangle of Figure 13.11 we have ^/rj/ yfr = 
sin(7r^)/ sin(7r r ). □ 



3.2 Growth model in the anisotropic KPZ class 

For fixed r, the macroscopic slopes of the interface in the v- and ^-directions 
are given by h v := d u h and h v := d v h. The speed of growth of the surface, 
d T h, depends only on these two slopes. Indeed, by (13.151) . we can rewrite 



dh 1 sm(Trh,,) smiirhr, 



(3.22) 



Or 7r sin(7r(/i v + h^)) 
Remark that the speed of growth is monotonically decreasing with the slope 



dh v 

for h v , h v , h v + h v e (0, 1). 



dv(h u , hv) <Q dv(h u , h^) ^ Q 



dh n 



(3.23) 
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To see which universality class our model belongs to, we need to compute 
the determinant of the Hessian of v — v(h l/ , h^). Explicit computations give 



d K d K v d K d hn v 



2 M^) 2 M^) 2 < (3 24) 



for hy^h^.hy + hrj e (0, 1), i.e., for (is, r], r) 6 P. Thus, our model belongs to 
the anisotropic KPZ universality class of growth models in 2 + 1 dimensions. 

3.3 A few other geometric properties 

During the asymptotic analysis we will use a few more geometric quantities, 
which we collect in this section. The key function to be analyzed is 

G(w) = G(w\v, rj, t) — rw + i/ln(l — w) — r]ln(w), w6C. (3.25) 

The critical points of G coincide with Q as stated below. 

Proposition 3.8. On C \ {0, 1} ; the function G has two critical points 
(counted with multiplicities). These two points are distinct and complex con- 
jugate if and only if (u, r], r) e T>, in which case the critical points are {Q, Q}. 

Proof of Proposition \3.8l The derivative of G gives 

n't \ T (( V + r - u\ 2 Ar]T - (n + t - uf\ 

GH = M^I) [{ W -—^—) + 4^ J' (3 ' 26) 

and we have two distinct complex conjugate solutions iff Atjt— (t]+t— v) 2 > 0, 
i.e., iff (u, rj, r) G V. Also, from (13.51) and (13.61) we get 



H. ( n) = 2±izj:, Ma) = y&EHn*?. ,3.27) 

It It 

Thus, Q and Q are the two solutions of G'(w) = 0, i.e., the two critical 
points. □ 

The main formulas needed later are the partial derivatives of Q as well 
as G"(ft). 



Proposition 3.9. Denote k = 2rlm(0) = a/4??t — (rj + r — v) 2 . Then we 
have 

G "< n > = «ray < 3 ' 28 > 
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which implies 



= 10(1-0)1 ' ^"{{1)) = -*--^ + ^. (3.29) 



Moreover, 

dn _ iO 90 _ i(l - O) dO _ -iO(l - O) 

dv k : dr] k ' dr k 

Proof of Proposition \3.9l From (13.26?) we get 



(3.30) 



The modulus is immediate, while the argument is obtained using (13. 5p . 

Since O is the intersection point of the circles \z\ = \[t\/t and 
|1 — z\ = v/ v/t, the direction of d v Q. is orthogonal to the vector and d v il 
is orthogonal to 1 — 0. Therefore, for some C\, c% G R, 

dfl <90 

— = Cl 0i, — = c 2 (l-0)i. (3.32) 

Looking at the real part of these equations, we get ctRe(O) = — cilm(O), 
and <9t,Re(0) = C2im(0). On the other hand, 

Re(O) = - + o T - - ^Re(O) = -— , <9„Re(0) = — . (3.33) 

2r 2r 2r 



From this we conclude that 



2rlm(0) ' 77 2rlm(0) 



d » Q = o-t_/^ » = o_t^/ov ( 3 - 34 ) 



To get <9 T 0, we can use the following property: 0(av, ar/, ar) = 0(z/, 77, r) for 
any a > 0, which implies 

+ ^ + r<9 T ) = 0. (3.35) 

This equation leads to 

W = F~~7rvT f-n + ^(1 - O)) = - iQ( T 1 ~" ) , (3.36) 

2rlm(0) Vr r v V 2rlm(0) ' v ; 

using |0| 2 = rr/r and |1 — 0| 2 = v/t, see (13. 51) . □ 

Another important function appearing in the asymptotics of the kernel is 
the imaginary part of G(O) (and its derivatives). 
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Proposition 3.10. We have 

7(1/, 77, t) := lm(G(Q)) = rlm(fi) — vix r] — r]n u , (3.37) 
Its derivatives are 

Wv = _7r " d~r, = "^' (3 ' 38) 

and 

d 2 Im(G(tt)) 1 , n . , 

V = > k — 2rlm (fi . 3.39 

ouorj K 

Proof of Proposition ^ . ICk The relation (13.371) is a direct consequence of (13.51) . 

The rest are just simple computations. □ 



4 Gaussian fluctuations 

In this section we first state a couple of equivalent forms of the correlation 
kernel. In particular, the kernel for a fixed (n, t) has a Christoffel-Darboux 
representation in terms of Charlier polynomials. In the second part of the 
section we prove Theorem 11.21 on the Gaussian fluctuations. 

4.1 Kernel representations 

For the analysis of the variance we will use a representation in terms of 
Charlier polynomials. These polynomials are defined on Z + = {0, 1,2,.. .}, 
while our particles at level n live on {— n, —n + 1, . . .}. Thus, it is convenient 
to shift the position at level n by n, i.e., the positions of particles at level 
n will be denoted by — n +x, x > 0. Finally, we also conjugate by a factor 
(— l) ni_n2 . More precisely, the relation between the shifted and conjugate 
kernel K and the kernel /C in Theorem 11.11 is the following, 

K(x u m, ti; x 2 , n 2 , t 2 ) = (-l) ni ~ n2 }C( y x 1 - n u n u ti, x 2 - n 2 , n 2 , t 2 ). (4.1) 

For later use, we give the explicit double integral representation of K 
which will be used in the asymptotic analysis. 

Corollary 4.1. The extended kernel K can be expressed as 
K(x 1 ,n 1 ,t 1 ;x 2 ,n 2 ,t 2 ) 

[ w &i dz &o dw e^ ( i-V +1 et2W[ X w)X2 ^. K *i) ^ K t 2 ) 

( I2W dz 3>r . dw £tiz(1 _, ri+1 , K, *i) ~< K t 2 ) 

(4.2) 
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Proof of Corollary \4-l\ The kernel (14.21) is obtained by substituting into 
(14.11) the expression for K, from (11.61) . and applying the change of variables 
z ->• 1/(1 - w) and w -»• 1/(1 - z). □ 

It is instructive to see the structure of the kernel that leads the above 
integral representation. 

Proposition 4.2. The extended kernel K is given by 



k=l 



(4.3) 



with 



^ n '\x) = — <fi dw - 

,n,t, v -If, z x e~ t: 



i 



e tw (l-w) k 



•» <*> = 2iif ri d! (M^' (44) 

i r pw(t!-t 2 ) 

0((ni,ti),(n 2 ,t 2 ))/ s = A d IrfnitiWfnatall, 

r v zy 27Ti Jp tyai-aa+l (^ — l)n 2 -ni U n lW^tn 2 ,r 2 ,|j ? 

where T^i and Ti are any simple anticlockwise oriented contours that include 
poles {0, 1} and {1}, respectively. 

Proof of Proposition \4-2\ Using the integral representations for \1/ and $ one 
checks that 

^2^l lM (x)^ t2 (y) = ^ ni ^ n ^\x,y). (4.5) 
Thus ( 14. 3 p becomes 

E?=iK\%(xi)Kl'%(x2), (n utl ) * (n 2 ,t 2 ) 

(4.6) 

This new expression is good because in (14.4ft we never have the case when 
the pole at w = 1 in ty^' survives. Then, one has just to take the sums inside 
the integral. For example, for (ni,ti) (n 2 ,t 2 ), we first take the sum inside 
the integrals and then we extend it to k = oo. This can be done provided 
|1 — w\ > |1 — z\. Then, to get the formula (14. 2p . one just has to rename the 
variables z — > 1 — w and w — > 1 — z. □ 

For the computation of the variance, we will need only the kernel at fixed 
(n,t). It is given in terms of the Charlier polynomials, C n (x,t), given by 

c „ (M) = ^/ dm (i^, (4 . 7) 

t n 2711 / r „ W n+i 



K(x 1 ,n 1 ,t 1 ;x 2 ,n 2 ,t 2 ) 



ii 
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which satisfy C n (x,t) = C x (n,t), and are orthogonal with respect to the 
weight Wt(x) = — p, namely 

2^ C n (x, t)C m (x, t)w t (x) = —S n ,m, t > 0. (4.8) 

Corollary 4.3. The kernel K(x,n,t;y,n,t) is equivalent (conjugate) to the 
kernel K n j(x,y) , given by 

is t \ r-r ( ln-i(x,t)q n (y,t) - q n (x,t)q n _ 1 (y,t) 

K nt {x,y) = Vnt , (4.9) 

x-y 

where 

i t n/2 

q n {x, t) = Wt ( x fl 2 C n {x, t), (4.10) 



Proof of Corollary \4.3\ Consider n\ = n% = n and t\ = t 2 = t m. (14.31) . Then, 

n-1 

K(x, n, t; y, n, t) = £ ^'(x)$f (y). (4.11) 



fc=0 

For all k > 0, w = 1 is not a pole in the integral representation of ^' t - 



Using g2D and C n (x,t) = C x (n,t), we get ^\x) = %C k (x,t). Also, by the 
change of variable z — > 1 — w in the integral representation we obtain 
^fc'*! 3 ') = e < ^Cfc( a; ) ^)- Thus the kernel writes 

K(x,n,t;y,n,t) = w t (x) Y,jfk{x,t)C k {y,t), (4.12) 

fc=0 

which is conjugate to the kernel 



n— 1 



#n,t(z, 2/) = XI qk ( x > ^(y, t). (4.13) 

fc=0 

From C n (x,t) = u n x n + - ■ ■ with u n = l/(— t)™, we have q n (x,t) = v n x n + ■ ■ • 
with v n = (—l) n /^/t n n\. Then, (14. 9 j) follows from Christoffel-Darboux for- 
mula. □ 

Remark 4.4. For later use, we rewrite q n as 

e -t/2 t x/2 

g„(x,t) = B ntt (x)I ni t(x), B ntt (x) = tn/2 , (4.14) 



and 

I n J x ) = —l dw ^~ W iy W \ (4.15) 
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4.2 Proof of Theorem Q 

In this section we look only at the height function at a given time. Therefore, 
it is convenient to set A = vjr and c = rj/r so that we have n = [rjL] = [ct] 
and x = [uL] = [At]. In these variables, the equation for the bulk region 
given by (13.31) rewrites as 

(1 - v^) 2 < A < (1 + v^) 2 . (4.16) 

First we compute the variance of the height. 

Proposition 4.5. For any A G ((1 — \fc) 2 , (v^+ 1) 2 )> 

Var(fr([(A-c)t],[ct],t)) = J_ 
i-oo ln(t) 2tt 2 ' 1 ' J 

With this we can prove Theorem II .2[ 

Proof of Theorem {TM It is a consequence of Proposition 14.51 and [50]. More 
precisely, in Section 2 of [50] the convergence in distribution (a generalization 
of the result for the sine kernel of [18]) is stated. However, following the proof 
of the theorem, one realizes that it is done by controlling the cumulants, i.e., 
also the moments converge. □ 

Proof of Proposition "J75, The variance can be written in terms of the one 
and two point correlation functions p\ and pi- Namely, 



2 

[X) ' 



V&r(h([(\-c)t],[ct],t))= Yl P2(x,y)+ ^ Pi(x)- ( J] pi(a 

x,y>[\t] x>[\t] x>[Xt] 

(4.18) 

where p 2 (x,y) = K n>t (x,x)K n;t (y,y) - K n)t (x,y)K n>t (y,x) and px(x) = 
K n>t (x,x). Using K* <t = K n>t on 4(Z + ), we have 

Var(/i([(A - c)t], [ct],t)) = ^ K n ^(x,x) - £ K n , t (x,y)K nit (y,x) 

x>[\t] x,y>[\t] 

oo 

= JZ^*^'^)^ 71 '*^'^ ~ K nA x ^y) K n,t(y^ x ) 

x>[Xt] y=0 x,y>[\t] 

= ( K nA x >y)) 2 > n=[ct]. 

z>[At] i/<[At] 

(4.19) 
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We use the expression ( 14.91) for the kernel K n>t . We decompose the sum 
in (I4.19P into the following three sets: 

M = {x,y eZ 2 + \x > [Xt],y < [Xt],y - x < ext}, 

Ri = {x,yeZ+\x>[\t],y<[Xt],e 1 t<y-x<e 2 t}, (4.20) 

R 2 = {x,y e Z 2 + \x > [\t},y < [\t],e 2 t < y - x}, 

where the parameter £2 = 5 m i n {(l + v^) 2 — A, A — (1 — yfc) 2 } is chosen so 
that R\ is a subset of the bulk. Thus 

Var(/i([(A - c)t], [ct],t)) =M t + R t<1 + R^, (4.21) 

with 

M t = Y \K n , t (x,y)\ 2 , Rt,k= Y \Kn, t &y)\ 2 ■ ( 4 - 22 ) 

x,y£M x,y£R k 

Remark: The parameter E\, small, will be chosen i-dependent in the end. 

(1) Bound on R t)2 . For x, y G R 2 , we use y — x > e 2 t, and extend the sum 
to infinities 

Rt,2 < \ Y Y (hct]( x ,t)\ 2 hct]-i(y,t)\ 2 + \q[ct]-i{x,t)\ 2 \q [ct] {y,t)\ 2 \ 

2 x>\t y<\t 

4 

+2\q [ct y 1 ( y x,t)q [ct ](x,t)\\q [cthl (y,t)q [ct] {y,t)\) < — . (4.23) 
The last inequality follows from Cauchy-Schwarz and the property 

Y\Qk(x,t)\ 2 = W*, = forallfc - ( 4 - 24 ) 

x>0 

(2) Bound on Rt^. Since this time x,y G R\ are always in the bulk, we 
just use the bound of Lemma [6781 and get 

V- 1 ^ 1 

R t 1 < const > -. rr = const > - 

' - ^ (x-y) 2 ^ z 

x,yeRi v y/ z=[eit] 

= *([e 2 t] + l) -*([£!*]), (4.25) 

where ^>(x) is the digamma function, which has the Taylor expansion at 
infinity given by 

(x) = ln(x) - l/(2x) + 0{l/x 2 ). (4.26) 

Thus 

Rt,i < const ln(l/£i), (4.27) 
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with const t- independent, as long as z, t — > oo as t —>■ oo. 

(3) Limit value for M t . This time we need more than just a bound. Recall 
that n = ct and set x = [At] + £1, y = [At] — £ 2 - We have 1 < £1 + £ 2 < £it- 
Lemma 16.41 gives 



£ -l/2 



rr 4/^ _ (l+c-A) 2 



0(t" 1/2 ) + 0(ei 



+ cos 



ta(c, A + £/t) + P(c, A) - £d c a(c, A) 



(4.28) 



We use it with £ = 0, 1, together with the trigonometric identity 

cos(6i + 5) cos(6 2 ) — cos(fei) cos(6 2 + 5) = sxn(6) sin(fe 2 — (4.29) 

with 5 = — d c a(c, A), 61 = ta(c, A+£i/t)+/3(c, A), 6 2 = ta(c, A— £i/t)+/3(c, A). 
The factor sin 2 (5) cancels the term exactly. We obtain (using (14. 9p ) 



6=16=0 



vr 2 (6 + 6) 



+ sin 



t(a(c, A - £ 2 /t) - a{c, A + 



(4.30) 



The contribution of the error terms can be bounded by hi{e\t)0{t 1//2 ,£i) 
and the remainder is 



[eit] Ci-1 

EE 



1 1 

^(6 + 6) 2 



sm 



t(a(c, A - £ 2 /t) - a(c, A + . (4.31) 



6=1 6=0 
Let 6(A) = — a(c, A), then 



b'(X) = arccos ° ) e (0, vr), for (1-y^) 2 < A < (1 + v^) 2 - (4.32) 

By Lemma 14.61 below (one needs to shift the argument of b(X ± £j/t) by A 
to apply it), for large t the leading term in the sum is identical to the one 
where sin 2 (- ■ • ) is replaced by its mean, i.e., 1/2. Thus 



sit] 6-1 



(H3U) 



6=16=0 VS1 s 
^ln^a + O^,^ 2 ^)- 1 )). 



(4.33) 
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Thus, 

M t = In( ei t) ( JL + ^r 1/2 , £i, (^V*)- 1 )) • (4.34) 
Now we choose E\ = l/ln(t). Then, 

Vax(h([(\-c)t],[ct],t)) = -^ln(t) + 0(l,ln(ln(t)),(ln(t)) 3 M), (4.35) 

which implies (14.171) . Modulo Lemma [4.61 the proof of Theorem 11.21 is com- 
plete. □ 

Lemma 4.6. Letb(x) be a smooth function (C 2 is enough) on a neighborhood 
of the origin with fe'(O) G (0, 7r). Then 

^1 ^' 3 in 2 [afc/t) - M-hIM _ y y 1 i (, , _J_ 

(4.36) 

uniformly for e > small enough. 

Proof of Lemma \4-6] We divide the sum into two regions 

h = {&> 1,6 >0|l<& + 6 <£</*}, (4.37) 
/ 2 = {6 > 1, 6 > 0\eV~t < 6 + 6 < 

Let us evaluate the contribution to (14.361) of (£i> 6) £ ii. We set z = £1 + £2 
and get 

E E -2 sin2 - - *)/*)] • ( 4 - 38 ) 

2=1 g! = l 

Taylor expansion around zero leads to 

tb(^/t) - - «)/*) = z6'(0) + 0(£ 2 ). (4.39) 

Thus 

g3SD = ^ i (sin 2 [z6'(0)] + £>( £ 2 )) . (4.40) 

z 

The sum with the sine squared can be explicitly evaluated: 

^ sin>z) = l ]n(p) + 0(1)| as p^oo (4.41) 

2 = 1 
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provided < a < vr. Since J2z=i l l z = H p )/ 2 + 0(1/ P), we have 

-^T^ = ^ - (1 + 0(1/P)) . (4.42) 

2 = 1 " 2 = 1 

Using P = [e\/t] and going back to the original variables (£1,^2) we have 



sm 2 [tb(Zi/t)-tb(-&/t)] ^ 1 ^ , 1 £2 



(4.43) 

Now we evaluate the contribution to (I4.36P of (^1,^2) £ h- Let 
(X, Y) E I2, then we have X + Y > e\/t. We consider a neighborhood of size 
M = [e 2 \/t] around (X, Y), namely the contribution 

M 

E (x + y + a . + v) 2 sin2 W x + *)/*) - *(-(^ + y)/*)] • ( 4 - 44 ) 

a;,i/=0 

Since sin 2 (- ■ ■ ) > and - {x+Y \ x+y) 2 > 0, if we replace {x+Y l +x+y) 2 

by ( X +Y) 2 m (I4-44P the error made is bounded by 

M ' 1 1 



E(( X + y)a (X + Y + x + yy) (445) 
E fy±vi2 ( 1 ~ n Z rv7vi2 ) = E 1 vW ^)' 



„(*-nH (i + o(.)) 2 ; 

because (x + y)/ (X + Y) < 2e. This relation can be inverted and we also get 
m 1 m 1 

E(YW^ (A- + K + . + y p ' 1 + O ^)' < 4 ' 46 ) 

x,j/=o %,y=o 

Therefore we have 

ra = V °( e ) i V sin 2 [t&((X + x)/t)-t&(-(r + y)/t)] 

^ n (x + r + x + y) 2+ ^ n (x + Yy 

x,y=0 v ' a;,y=0 v 

(4.47) 

Now we apply Taylor expansion to the argument in the sine squared. 
Denote k x = t&(J¥/t) - tb(-Y/t), 6 1 = b'(X/t) and # 2 = b'(-Y/t). Then the 
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argument in the sin 2 (- ■ ■ ) is K\ + Q\X + 6 2 y + 0{e 2 ). The e 2 error term is 
smaller than the 0(e) in f!4.47p . thus 

dOTD=V ^ Asin^H^ 

x,y=0 v y x,y=0 

Since 6 is smooth and b'(0) G (0,7r), in a neighborhood of we also have 
b' G (0,7r). Thus, for e small enough, < 9i,9 2 < n uniformly in t, because 
\Y\/t < e and \X\/t < e. The second sum in (14.481) can be computed 
explicitly. For < 9i, 9 2 < it we have the identity 

M 

2j sin 2 [ki + Q\x + 9 2 y] 

x,y=0 

(M + l) 2 _ cos(2/ tl + 9iM + 9 2 M) sin(fli(M + 1)) sin(fl 2 (M + 1)) 
2 2sin(#i)sin(# 2 ) 

= E 2 (1 + ° (1/m2)) - (4 - 49) 

x,y=0 

We replace f !4.49j) into (14.48j) and finally obtain 

A sin 2 [tfe((X + z)/t) - tb(-(Y + y)/t)} 
^ (X + Y + x + y) 2 

x,y=0 v ' 

M 

= E 2 ^ + y + J + y) , (i + ^.^)-))- 



(4.50) 



This estimate holds for all the region ^2, thus 

V- sin 2 [^(6/t)-t&(-6A)] _ 1 ( , 1 ^ g r^-^ 

(41,42)6/2 (f 1,52)6/2 

(4.51) 

The estimates of (14.431) and ( 14.511) imply the statement of the Lemma. □ 



5 Correlations along space-like paths 

In this section we present an extension of Theorem 11.11 to the three types of 
lozenges (see Figure l5TTj) . Then we explain the three different ways of com- 
puting height differences. These are then used in the proof of Theorem 11.31 
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Type of lozenges 



Facet's types 




O = (x, n) position 



Type / related angle 



I / ix r] 



H/vr r 



III / 7T, 



v 



Figure 5.1: Facet's types of Figure fl~2l their associated lozenges and angles. 

5.1 Joint distribution of the three types of lozenges 

As we saw in the introduction, particle configurations can also be interpreted 
as lozenge tiling (see Figure ll.2p of a half-plane. One can draw the corre- 
sponding triangular lattice by associating to the three types of facets three 
lozenges made by one black and one white triangle as indicated in Figure 15.11 
We define the position of a black / white triangle to be (x, n)-coordinate on 
the mid-point of its horizontal side. 



Thus, in our system of coordinates, these positions are pairs of integers (x, n) 
with x G Z, n e {0, 1, . . .} for black and n 6 {1,2,.. .} for white triangles. 
We first state the result in the common way from tiling point of view, and 
then we will reformulate it by using the kernel K defined in ( 14.21) . 

For any pair of black and white triangles with space-time coordinates 
(x,n,t) and (x',n',t'), define the kernel 



/C(L (x, n,t);\\ (x', n', t')) = (_i)*-* , +»-n' K.{x, n, t; x' , n, if), (5.1) 



where K, is the kernel defined in (II. 6ft . 

Theorem 5.1. Consider a finite set of lozenges at time moments t\ < ti < 
■ ■ ■ < 1m, consisting of triangles 




(bi,Wi) := (^(2r<,ni,*i),\|(a4nJ,fi)). 
Assume that > nj if ti < tj. Then 

P^There is a lozenge (bi,Wi) at time ti, 

for every i = 1, . . . , M} = det /C(6«, Wj) 



-I 1<M<A/ 



(5.2) 



(5.3) 
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Proof of Theorem I5.il We prove the statement by induction on the number 
of lozenges (5j, Wj) which are not of the form ^ . When this number is zero, 
then the statement is Theorem II .![ which is the base of the induction. 

Consider any set S of lozenges at any time moments, plus another lozenge. 
Then, the l.h.s. of ( 15.31) obviously satisfies 

P(5u^)+P(5u\k)+P(5UB) = P(5), (5.4) 

where in the l.h.s we either keep the white triangle fixed, or we keep the black 
triangle fixed (and we assume that S does not contain the fixed triangle). 
Next we verify that the same relation holds for the r.h.s. of (15.31) . 

Case (a): the fixed triangle is white. From the explicit formula for the 
kernel, we get 

K(k {x,n,t);X\)+jC(± (x,n-l,t);^)+)C(± (x + l,n - l,t);X\) 

= f 1, if^ =S\(x,n,t), 
\ 0, otherwise. 

(5.5) 

This implies relation (15.41) for the r.h.s. of (15. 3p . 

Case (b): the fixed triangle is black. There are two possibilities: (i) 
the black triangle is on the lower boundary {(x, n,t)\n = 0}, or (ii) it is not 
on the boundary. In case (ii), the relevant relation on the kernel is 

£(k ; ^ n',t))+IC(^;X\ (x', n ' + l,t)) + JC{^ ; ^ {x 1 - 1, n' + 1, t)) 

= f 1, ifk = k.(x',n',t), 

\ 0, otherwise. 

(5.6) 

In case (i), our assumption rii > n'- whenever ti < t'j implies that we are con- 
sidering the last time moment, tu- Then, in the formula for the kernel (11.61) . 
the first residue term drops out and the second term vanishes on \| (a/, 0, t) 
(since at z — 1 there is no pole anymore). Thus ( 15. 6ft still gives the needed 
relation: 

P(5U^k)+P(5uB) = P(5). (5.7) 

With the relation ( 15. Ah verified (which, in one case, degenerates to 
( 15.71) ). let us explain the induction step. Let us take a lozenge in the set 

{(bi,Wi),l < i < M} which is not of the type ^ . For example, consider 
B and denote by S the set of remaining M — 1 lozenges. Then 

P(5UB) = P(5)-P(5U^)-P(5U^), with ^ fixed. (5.8) 
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So, we have a linear combination of two terms with one less lozenge of type 
different from ^ , plus the third term with XL. whose black triangle is one 
position on the right with respect to the B . For this term we use 

P (S U XL) = P (<S)-P (s U ' ^) -P (S U \E) , with ^ fixed. (5.9) 

So, the third term in (15.81) is rewritten as a linear combination of two terms 
with one less lozenge of type different from ^ , plus a third term with a 
lozenge of type B , and this lozenge is one position to the right from the 
initial B in the l.h.s. of ( 15. 8ft . This can be continued iteratively. A similar 
argument holds for lozenges of type XB with (15.81) and (I5.9P applied in the 
opposite order. 

Thus, we can represent the r.h.s. of (I5.3P as linear combination of those 
with fewer lozenges of type XB , B , plus an expression of the same kind but 
with one of the XBor B lozenges far to the right. 

We still have to verify that the formula with one more lozenge of type 
XBor B agree when such a lozenge moves to +oo. Since the determi- 
nant in (15. 3p is invariant with respect to conjugation, consider the kernel 
2 X_X '/C(B. (x, n,t);\\ (x', n', t')) instead. Then, one verifies that 

2 x - x 'jC(^(x,n,t);Xl(x',n,t')) -> as x' -»• +oo, 

(5.10) 

2 x- x ^ n ,t);X] (x', ri, 0) ^ as x ^ +oo, 
and also for a lozenge (b,w) far to the right (i.e., when x — ► +oo) we have 

Therefore, if in r.h.s. of (15.3p there is one lozenge XBthat is far to the right, 
the determinant tends to zero, which agrees with 

P(5UXW)^0 asXk^+oo. (5.12) 

On the other hand, if in r.h.s. of (15. 3p there is one lozenge Bthat is far to the 
right, the determinant tends to the determinant of its minor corresponding 
to S, which is in agreement with 

P(5UB)^P(5) asB^+oo. (5.13) 

This completes the induction step. □ 

In the next section we will describe height function differences as sum 
over lozenges of type ^ or XB . To each lozenge one can associate a position. 



60 



We decided to set the position of a lozenge to be equal to the position of 
the white triangle, see Figure |5~T1 Now we restate Theorem 15 .11 in a slightly 
different form. 



Theorem 5.2. For any N = 1,2, ... , pick N triples 

Kj = (xj, rij, tj) 6Zx Z>o x R>o 
such that Xj + rij > and 

ti < t 2 < • • • < t N , ni > n 2 > • • • > n^r. (5-14) 

T/ien 

P{For eac/i j = 1, . . . , N at (xj,Uj, tj) there is a lozenge 

of type e 3 e {I, II, III}} = det [Kefa, 6 t ; Xj, 0/)^ ,. (5.15) 

where 

Koix^nx, ti, Ox] x 2 , n 2 , t 2 , 9 2 ) 




K(xi + ni,ni,ti; x 2 + n 2 ,n 2 ,t 2 ), if 9 1 = I, 

— K(x\ +ni,ni — l,tx,x 2 +n 2 ,n 2 ,t 2 ), if Q\ =11, (5.16) 

K(x 1 +ni — l,ni - l,ti;x 2 + n 2 ,n 2 ,t 2 ), if 6 X = III, 



u^/j as defined in Section 

Proof of Theorem 1 5. M The proof is simple. One just applies the correspon- 
dence 

Type I at (x, n, t) (^ (x, n, £), \| (x, n, t)) (5-17) 
Type II at (x, n, t) <£> (^ (x + 1, n - 1, t), ^ (x, n, t)) (5.18) 
Type III at (x, n, t) & (^ (x, n - 1, t), X\ (x, n, t)) (5.19) 

to (15.31) and then rewrites /C in terms of /C. Then, using the relation (14.11) . 
we get the expression in terms of K. Finally, one conjugates the kernel by 
(— l) Xl ~ X2 and obtains the desired kernel Kg. □ 

5.2 Height differences as time integration of fluxes 

To determine the height function at a position (m, n) at a given time t, one 
can act in three different ways: 
(a) Sum along the x-direction: 

h(m, n, t) = 1 (lozenge of type I at (x, n, t)) . (5.20) 
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t 

p 

/ ) p -p 




Figure 5.2: The black dots represents the space-time positions where we 
want to study the height functions. They live on a space-like surface (i.e., 
for any two points (xi, ni, ti), (£2,^2,^2) on it, either (ni,ti) -< (712,^2) or 
(n2,t 2 ) -< (ni,ti)). The white dots represents the projection of the black 
dots onto the (n, t)-plane. 



(b) Sum along the n-direction: for n! > n, 

h(m, n, t) = h(m, n , t) + H njn r(m, t), (5-21) 

where 

n> 

H n)n i{m,t) = — 1 (lozenge of type II at (m,p,t)) . (5.22) 

p=n+l 

(c) Integrate the current over time: for t > t', 

h(m, n, t) = h(m, n, t') + J t > t t(m, n), (5.23) 

where Jt' it (m,n) is the number of particles (= lozenges of type I) which 
jumped from site (m,n) to site (m + l,n) during the time interval [t',t]. 

In principle, one could use (a) alone to determine the height. However, 
this turns out to be not very practical when dealing with joint distributions 
of height functions at different points (mi,n, t), . . . , (m^,n, t). The reason 
is that the height functions are linear functions of lozenges of type I but 
the same lozenges appears in several of them. The result is a very tedious 
computation. This can be avoided by using (b) and (c) depending on the 
cases, see Figure 15.21 for an illustration. 

Therefore, the expression 

E I Y[ [h(m k} n k , t k ) - E(h(m k , n k , t k ))} J (5.24) 
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can be expressed as a sum of terms of the form 

M 



k=l 



E( Y[[h(m k ,n k ,t k ) - E,(h(m k ,n k ,t k ))\ 

R 

IJ [H nttn r t (m e ,ti) - E(H nttn r t (mt,tt))] 



e=M+i 



j=R+l 



(5.25) 



We now derive a formula for (I5.25|) . 



Lemma 5.3. Assume that the following paths do not intersect and lie on a 
space-like surface: 



{(x,n k ,t k )\x > m k }, 
{{mi,p, t e )\p = n £ + 1, 

{(iTLj, Hj,t)\t G [t'ptj]}, 



n! t }, 



k = l, 
£ = M+l 
j = R+l, 



,M, 



Then 



, R, 
N. 



(5.26) 



xi>mi 



E E 

x M >m M p M+1 =n M+1 +l 



tR+l 



ds 



R+l 



R+l 



E 

PR=n.R+l 

A 1A A 1>2 A 1;3 

^2,1 ^2,2 ^2,3 

^3,1 ^3,2 ^3,3 



ds at det 



with the matrix blocks Aij as follows: 



(5.27) 



Ai t i — [(1 — 5ij)K(xi + Hi, n,i, ti] Xj + nj, nj, tj)\i<ij<M ' 
^2,1 — [Kyrrii + Pi,Pi — 1, Xj + rij, n 3 -, ^jJIm+ixi^-R, i<j<M ' 

■^3,1 \K{rfli -\- Hi, 7li, Si, Xj 4~ Tlj, Tlj, ^ j ) ] l< i < A^, Kj<M ' 
^1,2 — [-^(^i + Tli, n iy ti] 17lj + Pj,Pj, tj)\i<i<M, M+l<j<R ' 

A 2 , 2 = [(1 - 5ij)K(mi+pi,pi - hUimj +p j ,p j ,t j )] M+1 < iJ < R , 
A'3,2 — [Kyrrii + ^> n i, sf, rrij + Pj,Pj, tj)]ft + i<i<N, M+x<j<R > 
^1,3 = [-^A"(a?j + n i5 n i? t i5 m,- + n J; n i? s,-)] i< i < A f,i ? +i<i<Jv ' 
^2,3 = [<9s^K + p i; pi - 1, Ui m 3 + nj, n h Sj)\ M+1 < i < i i ijR+1 < i <^ • 
^3,3 = [-(1 - 6i,j)d s .K(mi + m, rii, sf, rrij + n jy nj , sj)] R+1 < i:j < N ■ 



(5.2* 
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Proof of Lemma \5.3[ Below we prove that 



, M R N x 

E( Y[h(m k ,n k ,tk) Yl H netn/e (m e ,t e ) JJ J t ^ t . (raj, rij) J (5.29) 
^fc=i e=M+i j=R+i ' 

is equal to (15.271) but without the 1 — 5ij terms in Ax ) i,A 2) 2, and A 33 . The 
fact that the subtraction of the averages is given by putting zeros on the 
diagonal is a simple but important property, which was noticed for example 
in [33] (see the proof of Theorem 7.2). 

For N = R, (I5.29P is a direct application of Theorem 15.21 to the formulas 
( I5.20p and fl5.22p . The absence of the minus sign in A 2) * is a consequence 
of the minus in the definition of the if's in (15.221) . Next we extend the 
result when N > R, by first considering N = R + 1 for clarity. Denote by 
7](x, n, t, 8) the random variable that is equal to 1 if there is a type 9 lozenge 
at (x, n, t) and otherwise. Recall that lozenges of type I are exactly what 
we call particles. The flux of particles can be written as 



D 

J t , t (m,n) = lim V\(m,n,Ti_i,I)(l - rj(m,n,Ti,T)) (5.30) 



with Ti = f + iAr, i = 0,...,D, At = (t-t')/D. The quantity 
rj(m, n, Tj_i, I)(l — r](m, n, Ti, I)) equals 1 iff the site (m, n) was occupied by a 
particle at time Tj_i and empty at time Ti. Each particle tries to jump inde- 
pendently with an exponentially waiting time. Every time a particle moves, 
it can also push other particles, but no more than one on each (higher) level 
n = const. In any case, since on each level there is a finite number of parti- 
cles, the probability that a particle has more than one jump during time At 
is (9 (At 2 ). Thus, the limit At — » is straightforward. 

To obtain ( 15.291) we have to determine expression at first order in At of 

E ^(m,n,r i _i,I)(l - 77(771,71,7*, I)) J^fKm^n^tj-, 0,-)^ • (5.31) 

Then, in the At — > limit we will get an integral from i! to t. 

Set K x , n (ti;h) = K' h ( x + n)^ t2 {x + n). Remark that in (lOi 

0(( n > T i)'( n ' r i-i)) (x, x) = 1. Then, since 7* > 7i_i, from (14. 3B we obtain 



(153TD = det 



K mj n(Ti-i,Ti-i) K m jJji-x\Ti) K e (m,n,Ti-i,I;q) 
1 - K m>n {Ti,Ti- X ) 1 — K mtn (Ti;n) -K e (m,n,Ti,l;q) 
K e {q]m,n,Ti- U l) K e (q;m,n,Ti,l) K(q,q) 

(5.32) 
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where by q we denoted the quadruples (rrij, rij, tj, 9j), for j G {1, . . . , Q}. The 
kernel Kg is a simple function of the kernel K, see (15.1 61) . The second line is 
just one in the diagonal minus the entries of the kernel (c.f. complementation 
principle in the Appendix of [14]). Written in terms of K it becomes as 
above, since the (2, 1) entry has a 1 coming from 0. Next we perform two 
operations keeping the determinant invariant: 



Second row — 
Second column — 

We get that (15. 32ft is equal to 



Second row + First row 
Second column — First column. 



det 



{fi-li Tj_i) ATd 2 K mtn (Ti_i; Ti-i) + 0(At 2 ) 

1 - 0{Ar) 0(At 2 ) 
K e (q; m, n, T;_i, I) Ard 2 Kg(q; m, n, 7i_i, I) + 0(Ar 2 

d 2 K m:n (T i _ 1 ; 7i_i) K e (m, n, r^i, I; q) 
d 2 K e (q;m,n,Ti-i,T) Kg(q,q) 



K e {m,n,Ti^l\ q) 
0{At) 
Ke(q,q) 



-At det 



£>(Ar 2 



where d 2 means the derivative with respect to r^-i in the second entry of the 
kernel. This formula and (15.301) imply 

/ Q 
E I J t/:t (m,n)Y\v( m j, n j,tj;O j ) 



ds det 



—d 2 Kg(m, n, s, I; m, n, s, I) Kg(m, n, s, I; rrij, rij, tj, 9j) 
-d 2 K e (m h rii, t h ^; m, n, s, I) K B (m h n iy t h 9^ rrij, rij, tj, 9j) 



The case of several factors J is obtained by induction. Expressing Kg for 
9 G {I, II} in terms of K only, and considering the fact that H is minus the 
sum of lozenges of type II, we obtain the result. □ 



(5.33) 



(5.34) 



5.3 Proof of Theorem E3 

Consider the expectation 



N 



El Y[[H m k,n k ,t k ) -E(h(m k ,n k ,t k ))} j. 



k=l 



(5.35) 



Our goal is to determine its limit as L —>■ oo under the macroscopic scaling: 

Vk + 



t k = r k L, n k = [r) k L], m k = [(y k - rj k )L], with v k G ((^% - V^) 
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Figure 5.3: The image of paths of Figure E2] under the map fi. 



The expression ( 15.35)) is a linear combinations of expressions from 
Lemma 15.31 The r.h.s. of (I5.27P contains an N x N determinant; let us 
write it as the sum over permutations a G Sn of terms each of which is sgn o 
times the product of matrix elements (i, ai), i = 1, . . . , N. 

The contribution of all permutations with fixed points is zero (because 
the diagonal matrix elements are zeroes). All other permutations can be 
written as unions of several cycles of length £ > 2. The contributions of 
the permutations with only cycles of length 2 lead to the final result, i.e., to 
prove the theorem we first need to show that the sum of the contributions of 
permutations with cycles of length i > 3 vanishes in the L — > oo limit. 

Consider a cycle of length £ > 3 and use the indices !,...,£ for the 
corresponding points (rrii, rij,tj). Let us order them so that 

T]i > T)<2 > . . . > T)£, T\ < T2 < ■ ■ ■ < Tg, no double points, (5.36) 

i.e., {jij,Tj) -< (77i-i,Ti-i)- 

For an ^-cycle we need to take the product of the kernels (or their 
time derivatives depending on the case), and do the summation over 
%k > [(~Vk + Vk)L], or over pe G [[rjtL] + 1, [^^]], or the integration over 
[t/L, TiL] depending on whether in ( 15. 27ft we have a sum or an integral. 

We first collect all the factors related to a fixed index i. There are three 
possible cases: 

(a) The index i is related to a sum over [(—T)i + Vi)L, oo). Then we have to 
analyze 

^2 K ( x > [ViL\,TiL] x ai )K(x a -i;x, [rjiL], nL)\ (5.37) 

x>[uiL\ 

(b) The index i is related to a sum over [[rjiL] + 1, [^L]]. Then we have to 
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analyze: 



^2 K([fiL]+P-[Vi L ]>P-' i -iTiI'',X(T i ) 



p=[ m L]+l 



x K(x a -i- [viL] + p- [riiL],p,TiL); (5.38) 



(c) The index i is related to an integrated variable. We have in this case 



The second kernel has a shift by one in the second x-entry. This comes from 
the identity 



which immediately follows from ( 14.2ft . 

We analyze these three expression in the L — > oo limit using results of 
Section 16.31 First of all, since w c — z c remains bounded away from zero all 
along the integrals/sums, the bounds of Section |6~31 imply that the contri- 
butions of the error term 0{L~ l l*) in (16361 are of the same order, namely 
(9(L -1 / 8 ). Therefore we can get rid of them immediately and we will not 
write them explicitly in what follows. 

Case (a) We divide the sum in three parts for which we use 
Propositions I6.9ti6.13l 

Case (a/1) Sum in the interval 




(5.39) 



d s 'K(m, n, s; m , n , s') = K(m, n, s; m! + 1, n , s'), 



(5.40) 




(5.41) 



Then, by Propositions 16. 1 2116. 131 




xeh 
const 



const 



exp 




(r,L)V3 ) 



x terms in x„., x -i 

°%i (J - 



(5.42) 



- LV3 



x terms in x„., x -i. 



Therefore, as L — > oo, the contribution of this sum goes to zero. 
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Case (a/2) Sum in the interval 

h = {x G N, {y/n + ^) 2 L - L 2 ' 3 < x < (yfn + ^) 2 L - iL 1 / 3 }. (5.43) 
By Propositions I6.11H6.121 

^K(x, [r)iL),TiL] x ai )K(>€ a -i;x, [^L],TjL) 

E const 
x terms m x CTi , x CT -i (5.44) 

xeh L^J rjiTi + rji - x/L) 2 

const 

< -JT/6~ X terms 111 ^r 1 - 

Therefore, as L — > oo, this contribution is also infinitesimally small. 

In the following (Cases (a/3), (b), and (c)) we will assume that all the 
entries x/s of the kernel are in T> and apply Proposition 16.91 and its Corol- 
lary [HUTU Let us justify it. The variables corresponding to time integration 
(Case (c)) and sum over the p variables (Case (b)) in (I5.27P are always in V. 
Therefore, the only x/s which are not in T> correspond to Cases (a/1) and 
(a/2) above. From Propositions 16.9116. 121 the contributions in the Xj variable 
are of order 

Oil) 



<\Jvi T i -\(Vi + T i- v i) 2 



(5.45) 



if Xj G V. The sum in Case (a/3) is then bounded by 0(1) because the sum 
is over O(L) sites and the square- root singularity is integrable. Even simpler 
is Case (b) where we never come close to the singularity and the sum is over 
O(L) sites. Finally, in Case (c), the integration is over a time span O(L). 
Therefore, the contributions of the terms of Cases (a/3), (b), and (c) are 
(9(1), and for every sum reaching the edge we get a factor 0(L~ 1/6 ). Thus, 
in the following we need to determine the asymptotics of Cases (a/3), (b), 
and (c) in the case where all the entries Xj's are in T>. 

Case (a/3) Sum in the interval 

h = {x G N, [ Vi L] <x<(^F i + ^qtfL - L 2 / 3 }. (5.46) 
Define the functions 

A(u, T], r) = = (5.47) 

27r|G''(OM,T))|y^F 

and 

F(u, v,r)=L Im(G(n(i/, V , r))). (5.48) 
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Then, by Proposition 16.91 we have 

} u K(x, [rjiL], nL; x ai )K(x a -i;x, [rjiL], nL) 

A(x/L,r] i ,r i ) 



L 



u(a t ) - uj(i) cu(i) - 



+ _1 e -2iF{x/L, Vi , n ) (5.49) 

u)(Gi) - uj(i) u>(i) - tufo-r 1 ) " 



H — e 

gi/3i(i) e -i/3 2 («) 

H h O 



2\F(x/L,r)i,Ti) 



x terms in xu., x -1 



where we used the notation u;(z) = f2(i^, r]i, Tj) and O means the other 12 
terms obtained by replacing u(<7i) by tD(<7j) and/or cu^a^ 1 ) by ^(a^ 1 ). 

First we want to show that the terms with F in the exponential are 
irrelevant in the L — > oo limit. For that, we sum over positions 
around any vL in the bulk. Then, for < x < L 1 ^ 3 it holds 

F(u + x/L, 77, r) = Ly(v, 77, r) + x^ 7 (z/, 77, r) + C(L^ 1/3 ), (5.50) 

where 7(1/, 77, r) = L _1 F(z/, 77, r). All the other functions (A, /3x(z), /3 2 (i), 
and u;(i)) are smooth functions in z/j, i.e., over an interval L 1 ^ 3 vary only by 
~ L~ 2 I 3 . Therefore we have to compute an expression of the form 

N-1 

J- e 2i ^+*/^V(z, + x/L, 77, r) (5.51) 

x=0 

where is a smooth function given in term of A, $2(1), and Thus 

N-1 

( I53TD = 0(z/, 77, r)e 2ii7( ^' T) - ^ e ite + 0{L~ 1 ' 3 ) (5.52) 
with 6 = 2d v ^{y, 77, r). Then, for < b < n, we use 

1 ^ pibN _ 1 

In our case, b is strictly between and tt as soon as we are away from the 
facet. When we reach the lower facet, b — > 0. However, in the sum over I 3 we 
are at least at a distance l?l 3 from the facet, i.e., b > const L -1 / 6 . Therefore 

[ (E53D 1 < const /(6JV) < const L~ 1/6 . (5.54) 
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Since this holds uniformly in the domain J 3 , we have shown that the con- 
tribution of the terms where the exp(±2iF) is present is at worst of order 
£-1/6 Therefore the only no n- vanishing terms in ( 15. 49ft are 



A(x/L,r)i,Ti 



,i/8i(») 



p -i/3 2 («) 



lo[i) uo(i) — uoia- 







X 



terms in x ai , (5.55) 



All the functions appearing now are smooth and changing over distances 
x ~ L. Thus, defining x = uL, the sum becomes, up to an error of order 
e>(L~ 1/3 ), the integral 



{y/n+^/m) 12 



dvA(v, T)i,Ti 



-i/9i(i) 



u(i) 



U0\O: 



3 i/3i(0 



-i/3 2 (i) 



+ 



CJ <7. 



x terms in x av (5.56) 



The final step is a change of variable. For the term with oo(i), 
we set = u(i) = Q,(i/,r)i,Ti). Denote the new integration path by 
= {f2(z/, r/i, Tj), v : {^/ri+ y/rji) 2 — > ^}. The Jacobian was computed in 
Proposition 13. 9[ namely 



dv 



lujyi) 



K 



2mA e^We-^iW. 



(5.57) 



For the term with uj{i) we set z i = uo{i) = Q(u,i]i, Tj) and r*_ = T l + . Then 
(I5.56P becomes 



27Ti ^ El 



€i=± 



dzt 



1 



+ O 



x terms m. a il a i 1 . 

(5.5? 



The factor —1 comes from the orientation of r*., see Figure |5~3"1 

Case (b) We sum in the n-direction from [rjiL] + 1 to [r/-L] for some !][ > r)i. 
While doing this, we do not exit the domain D remaining in the bulk. 
Therefore, the computations are just a small variation of the sum over J3 
of Case (a). The minor difference comes from the —1 shift in p in the entries 
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of the first kernel. By changing the variable a = p/L, we then obtain 

p=[ m L]+l 



lim } K([uiL)+p-[rjiL),p-l,TiL;H^)K(H-i-[uiL)+p-[rjiL),p,Ti 

p=[ m L]+l 
4 



daA(ui -rji + a, a, t {j 



LOK%) U)[l) — W (T, 



-1\ 



+ 



u((Ji) - oj(i) u(i) - uj(a i l ) 



x terms m %„.. x-\ 



(5.59) 

with uj{i) = Q(vi — rii + a, a, r»). For the term with lu(i), we set zf = u(i) and 
denote the new integration path by T l + = {Qfa — rji + a, a, Tj), a : rj^ — > rji} 
(we set the orientation of the path as in Figure l5~3l) . By Proposition 13.91 we 
get 

Ml = 1 = 2ix\Ae^e-^u(i)-\ (5.60) 
oa k 

The change of variable for the term with oj(i) is similar. The result is the 
same formula as (I5.58P (of course, with the new r*.'s). 

Case (c) The last case is when we do an integration over a time interval. 
Similarly to Case (b), we do not have to deal with the edges, since, by 
assumption, we remain in the bulk of the system. We need to compute 



TiL 



dtK([uiL], [r)iL],t] x at )K(x a -i; [u t L] + 1, [rjiL],t) 



tL 



&TA{Pi,T)i,T) 



-i/3i(i) 



uj(ai) - uj{i) uj(i) - uj(a i l ) 



k(i-w(<)) 



+ 
+ 
+ 



ufa) - uj(i) 






e i/3 2 (i) 


ufa) -Q(i) 




gift© 


e -i/9s(i) 



-2\F(ui,rn,T) 



rrr(l -u(i))e 



j-(l - uu{i))e 2iF ^T) 



u((Ji) - uj{i) - u{a i v ) 



x terms in x CTi , x ct -i. 



(5.61) 

The only rapidly changing function is F, which, as for the sum, makes the 
contributions of the term with it vanishing small as L — > oo. We do the 
same change of variable as above, i.e., zf = u(i) = 0(^,7^,7"). Denote the 
new integration path by = {flfa, rji, r), r G [t^Tj]}. The Jacobian is 
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computed in Proposition 13.91 namely 

d»Q = -i^)(l-^)) = _ 2niAe iH%-^ {1 _ w(i)) . (5.62) 
or k 

Thus, we obtain again (15.581) . 

Thus, after summing / integrating all the I variables, we get the contri- 
bution of the £-cycle, namely 

v ' ei,...,££=± 1=1 si e f 1=1 * 1 

= g£ e n^/ rl ^'--'/ rl d -ri5^^- (5 - 63) 

v ; ei,...,e £ =± i=l 1/1 1/1 ^ 2=1 i * — i 

where we set <To := er^. By Lemma 7.3 in [33], which refers back to [17], 



E II yTy = °» for£>3. (5.64) 

<r=£-cycle in iSf i=l 1 1-1 

Therefore, the sum of (15.631) over all possible ^-cycles on the same set of 
indices gives zero for £ > 3. 

We have shown that we have a Gaussian type formula (sum over all cou- 
plings) for points macroscopically away. We still need to compute explicitly 



the covariance for such points. The covariance is obtained from (15. 63ft for 
1 = 2. We need now to consider the signature of the permutation, which for 
a 2-cycle is —1. We thus obtain a sum of 4 terms which can be put together 
into (see the end of Section 7 in [33] too) 

y /•Q(i/l,?7l,TL) r^t{V2,rH,vi) y 

77T^2 / dzi / dz2 l ^ ( 5 - 65 ) 



47T 2 \(fi(l/*,7fr,7i) - ^(^2 ) ^2,T 2 ))(fi(^l,r7l,Ti) - f2(z/ 2 , 7/2, T 2 ) ) 

5.4 Short and intermediate distance bounds 

Let us first prove the short distance bound (ll.lSp . 
Lemma 5.4. For any Xj G T> and any e > 0, we have 

E(H L (x 1 )---H L (x N )) = 0(L t ), L^oo. (5.66) 
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Proof of Lemma 



5j 



Theorem 11.21 implies, for any integer m > 1, 



E(tf L (x,-) 2m ) = 0(ln(L) m ). (5.67) 
By Chebyshev inequality, 

> Xln(L)) = 0(1/X 2m ), P(|ff L (x,-)| > F) = G(ln(L) m /F 2m ). 

(5.68) 

The final ingredient is that |ifi,(xj)| = O(L), since on level n = L we have 
only L particles. Therefore, for any K, we can bound 

\E(H L (x 1 )---H L (x N ))\ < P(\H l (k 1 )\<Y,... } \H l (x n )\<Y)Y n 

+ P(3j s.t. |# L (x,-)| >y)0(L) JV 
= 0(r w ) + C(L Ar ln(L) m /r 2m ). (5.69) 

Taking K = L e l 2 and m ^> 1 large enough, we obtain 

|E(F L (x!) ■ • • H L (x N ))\ = 0(L £ ), for any given e > 0. (5.70) 

□ 

Next we give an intermediate distance bound, extending the result of 
Theorem 11.31 

Lemma 5.5. Consider the setting as in Theorem 11.31 If the points fli's 
are not closer than L^ l ^ 8N \ then the difference between the expectation 
E(# L (*i) • • • H L (x N )) and r.h.s. of (FTfy is 0(L- 1 /^). 

Proof of Lemma 15.51 It is a small extension of Theorem 11.31 For N = 1 
the two expressions are identically equal to zero. So, consider N > 2. We 
have \Qi — Qj\ > L -1 / 16 , so that the estimate of ( 16 .561) can still be applied 
in the proof of Theorem 11.31 All the error terms collected are (9(L -1 / 6 ) (see 
times at most N factors of order 1/|^ - 0,| = 0(L 1 /^). This 
accounts into an error 0(L~ l l 2A ). Now, since the Qj's are not away of order 
one, when one 0(L~ 1 / 8 ) in ( 16. 56ft is used, it has to be multiplied by at 
worst N — 1 factors of order l/|Of — Qj\ = 0{L l ^ 8NS) ). Therefore the error 
is at most 0(L^ N - 1 ^^ L~ 1 ' 8 ) = O^L' 1 ^). Similarly, the contribution 
where (9(L~ 1//8 ) is used n times is 0(L^ N ~ n ^^ 8N ^ L~ n / 8 ), which is maximal 
at n = 1. Thus, for iV > 2, we get all together 0(L~ l l 2A ) + 0(L~ 1 /^) = 
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5.5 Gaussian Free Field 



The Gaussian Free Field on H, see e.g. [49], is a generalized Gaussian process 
(i.e. it is a probability measure on a suitable class of generalized functions on 
H) that can be characterized as follows. If we denote by GFF the random 
generalized function and take any sequence {4>k} of (compactly supported) 
test functions, the pairings {GFF(0k)} form a sequence of mean normal 
variables with covariance matrix 



E(GFF(0 fc )GFF(0j))= / \dz\ 2 {V<j> k (z), V<j>i{z)) 

Jm 

= / \dz 1 \ 2 \dz 2 \ 2 <f)k{zi)4>i{ z 2)G{zi,Z2) 
Jm 2 

where 



Q(z,w) = -^-ln 



z — w 



z — w 



(5.7i; 



is the Green function of the Laplacian on EI with Dirichlet boundary condi- 
tions. 

The value of GFF at a point cannot be defined. However, one can think 
of expectations of products of values of GFF at different points as being finite 
and equal to 

E[GFF(z x )---GFF(z m )] 

if m is odd, 

J2 Q(z a{1) , Z a{2 )) ■ ■ ■ 0(z<7(m-l)» 2<r(m)) if m is even - ^^ 2 ) 
pairings a 

The justification for the notation is the fact that for any finite number of 
test functions, 

„ m 

E(GFF(0i) • • • GFF(0 m )) = JJ |d^| 2 0fc(^)E [GFF(^) • • ■ GFF(z m )] . 

m k=i 

(5.73) 

The moments (15.731) uniquely determine the Gaussian Free Field. 

To state the convergence results, we consider any (smooth) space-like 
surface W C M 3 in the rounded part of the surface. Namely WcP, and for 
any two triples (z/j, r/i, X;) G U, i = 1, 2, r)\ < r] 2 implies T\ > r 2 . 

Clearly, the mapping Q restricted to U is a bijection. Consider any smooth 
parametrization u = (wi, W2) of U. Denote by flu the map from u to H, which 
is the composition of the map from u to (z/, rj, r) and Q. Then, for any smooth 
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compactly supported test function / on U, we define 

(f,H L ) := / duf(u)H L (u), 
Ju 

where Hl(u) is as in fll.lGp . Then 

(f,H L )= [ \dz\ 2 J(z)f(Q u \z))H L (Q u \z)), 
Ju 

where J(z) is the Jacobian of the change of variables z — *■ u by Q t 



(5.74) 



(5.75) 



Theorem 5.6. For any m G N, and any smooth compactly supported func- 
tions fi,...,f m on U, 



lim E 

L— ►oo 



,k=l 



II \dz k \ 2 f^(z k )E [GFF( Zl ) ■ ■ ■ GFF(z m )] 

(5.76) 



fc=i 



where ff{z) := J^)/^ 1 ^)). 



Remark 5.7. Since moment convergence to a (multidimensional) Gaus- 
sian implies convergence in distribution, Theorem 15.61 implies that the 
random vector ((fk, Hl))™ = i converges in distribution (and with all mo- 
ments) to the Gaussian vector with mean zero and covariance matrix 

ll/Hl^l 2 (v/ fc H (^),v/f(z))|U, i=lj ... jm . 

Proof of Theorem \5.b\ We have 



E 



\\{fk,H L ) 



,fc=i 



E" 



J] |dz fc | 2 /=(z fc )E [H L {Sl u \ Zl )) ■ ..H L (StfM)] ■ 

(5.77) 



fc=i 



Theorem 11.31 and Lemma 15.51 allow us to determine the last expected value 
as soon as \z{ — Zj\ are away at least of order 5 := L~ 1 ^ 8m '. Denote by 



W™ = {(z u . . . , z m ) e H m s.t. \ Zi - Zj\ < 5, 1 < i < j < m}. 
Then, as L — > oo, we have 



(5.78) 



J] Idz^f^E [H L {Q u \ Zl )) ■ ..H L {n u \z m ))] 

S~ fc=l 
„ m 

/ IT \dz k \ 2 ff(z k )E [GFF(zi) • • • GFF(z m )] + ©(ZT 



(5.79) 



l/(12m)\ 



75 



Then, since the logarithm is integrable around zero (in two but also in one 
dimension), the L — > oo limit is simply given by 

„ m 

hm([52HD= / Y\\dz k \ 2 ff{z h )B[G¥¥{ Zl )---G¥¥{z m )\. ( 5 - 8 °) 



We still need to control the contribution coming from H m \ H™. Using 
Lemma 15.41 this is bounded by 



„ m 

/ II \d*k\*ft(z k )E [H L (pJ-{zi)) ■ ■■H L {n u \z m ))] 



< const 5 L 



2 re 



(5.81) 

where const depends only the functions fi, ■ ■ ■ , f m - Since 5 2 = L 1 /( 4m ) 
and e > can be chosen smaller than l/(4m), in the L — > oo limit this 
contribution vanishes. □ 

We actually have a stronger result. Indeed the same formula holds also 
for smooth functions living on one-dimensional paths. Consider now any 
simple path 7 on U and denote by s a coordinate on 7. Denote by Q y the 
composition of the map from s to (u, r], r) and Q, and by 71 C H the image 
of 7 by lly. Then, we define 



</,tf L > 7 := / dsf(s)H L (s) (5.82) 

J-y 

and we get 

(f,H L ),= [ dzJ^f^iz^H^iz)), (5.83) 



13 



with J 7 (z) the Jacobian of the change of variables from z back to s by fi 7 1 . 

Theorem 5.8. For any m G N, consider any smooth functions f\, . . . , f m of 
compact support on 7. Then 



lim E 

L— ►oo 



.fc=i 



„ m 

/ IT \tek\fi(z k )E [GFF( Zl ) ■ ■ ■ GFF(z m )] 

Aff fe =i 

(5.84) 



where fi{z) := J^z)f{n?(z)). 

Proof of Theorem \5.8\ The strategy is the same as in the proof of Theorem l5.61 
The main difference is that the contribution at small distances will be of order 
8L e . However, this is fine, since we can choose e < l/(8m). □ 
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6 Asymptotics analysis 



In this section we do the asymptotic analysis of the functions g n 's at the 
(upper) edge and at the bulk. These are used to obtain the Gaussian fluc- 
tuations in Section HI Then, we do the asymptotic analysis of the extended 
kernel in the bulk and provide some bounds at the (upper) edge, needed to 
prove the Gaussian Free Field correlations in Section [51 



6.1 Asymptotics at the edge 

First we will determine the upper edge asymptotic of I n>t defined in (14.151) . for 
which we apply exactly the same strategy as in previous papers (Lemma 16.11 
and 16.21 are almost identical to the computations of Propositions 15 and 17 
in [10]). We first explain the strategy and then give the relevant details. 

Lemma 6.1 (Upper edge). Let n = ct and x = (1 + \fc) 2 t + st 1 ^ 3 , for any 
c > 0. Then, 

lim tVaj (.,.) j~ff = ~k 2 M(k 2 s), (6.1) 
t^oo e~v ct (l + v c ) 

uniformly for s in bounded sets, with n 2 = c 1//6 (l + ^fc)~ 2 / 3 , and 
k 2 = (1 + ^/c) 1 ^ 3 c~ 1 ^ 3 . Here Ai(-) is the classical Airy function. 

Proof of Lemma \6.1\ The strategy is the following. With the replacements 
n = ct and x = (1 + \fc)H + st 1 ^ 3 in (14.151) . we have an integral of the form 



2vri J r 



I dze tfo(z)+t 2 / 3 fi(z)+t 1 ^f 2 (z)+Mz) 



(6.2) 

r 



for some functions fk(z), k — 0, 1, 2, 3. The s-dependence is in f 2 {z). 

Step 1: Find a steep descent patlfl for the function fo(z), passing through 
the double critical point z c given by the condition f'§{z c ) = fo( z c) = 0. In 
particular, the steep descent path will be chosen so that close to the critical 
point the descent it the steepest. Then, uniformly for s in a bounded set, the 
contribution coming from the integration path away from a (^-neighborhood 
of z c is of order (9(e~ M *) with [i ~ 5 3 . 

Step 2: Consider the contribution of the integration only on \z — z c \ < 5, 



1 For an integral / = f dze^^', we say that 7 is a steep descent path if (1) Re(/(z)) 
reaches the maximum at some z £ 7: Re(/(z)) < Re(/(z )) for z e 7 \ {^o}, an d (2) 
Re(/(z)) is monotone along 7 except at its maximum point zq and, if 7 is closed, at a 
point z\ where the minimum of Re(/) is reached. 
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with 5 which can still be chosen small enough, but t-independent. In a neigh- 
borhood of the critical point, we can use Taylor expansion of the functions 
fo, . . . , / 3 and get 

eMtfo(zc) + t 2l *h{z c ) + t l ' 3 f 2 {z c ) + h{z c )) 

x -^r / dzexp(tK, (z - z c f/3 + t^K^z - z c ) 2 + t l,3 K 2 (z - z c )) 

2tti Jr n\z-z c \<5 

x exp(0(t(z - z c )\t 2 ' 3 {z - z c )\t l '\z - z c ) 2 , (z - z c ))). (6.3) 

Remark that we do not have a term t 2 ^ 3 (z — z c ) in the exponential. If such 
a term remains, than the edge scaling in x is not the right one. 

Step 3: Estimate the error terms. We do the change of variable 
t 1 / 3 ^ — z c ) = w and choose 5 small enough, so that the error terms are much 
smaller than the main ones. Subsequently, taking t large enough, the cubic 
term dominates all the others. Applying | e 2 ' — 1 1 < \y\e^ y]{ with y standing for 
the error term 0(- ■ ■ ), and changing the variable t l l 3 (z — z c ) = w, one sees 
that the difference between the integral with and without the error term is 
of order ©(r 1 / 3 ). 

Step 4: For the integral without errors, we also do the change of variable 
t 1 / 3 ^ — z c ) = w and then we extend the integration paths to infinity. This 
accounts for an error of order C(e~ M *). The final formula is then 

t-V3 exp(tf (z c ) + t 2 ' 3 h{z c ) + t^ 3 f 2 (z c ) + h{z c )) 

X / dwe K0U,3/3+Klw2+K2W + 0(r 1/3 ,e-^ , (6.4) 

where the integral goes from e _7ri//3 oo to e wl / 3 oo if k > and, in case kq < 
it goes from e~ 2m l 3 oo to e 27rl / 3 oo. The sign ±1 depends on the position of 
the critical point: we have +1 if z c > and —1 if z c < 0. Finally, the integral 
can we rewritten in terms of Airy functions using the following identity 

/ e az3/3+bzl + cz dz = a -l/3 Ai(fe 2 /fl 4/3 _ fi/ ^3^/3^ -bc/a ^ 



1 

2tx\ 



Specialization to our case. In our specific situation, the critical point is 
z c — —\/c, and the functions f , . . . , f 3 are 

fo(z) = z+{l + ^c) 2 \n{l-z)-c\n{z), 

fi(z) = 0, 

f 2 (z) = sln(l-z), 

f 3 (z) = -\n(z). (6.6) 



78 



Figure 6.1: The steep descent path used in the asymptotic analysis is the 
bold one. 



The steep descent path used in the analysis is made of pieces of the two 
following paths, 7 P and 7i oc (see Figure IBTTj) . given by 

7p = {-p^,<f> G ( 7r, 7r] }, 7ioc = {-V^ + e-^ 3B ^\x\,x G [0,^/2]}. 

(6.7) 

For p G (0, y/c\, 7 P is steep descent path for / . In fact, we get 

dRe{f {z = pe^)) psin0 2 , n r o a\ (a q\ 
d0 = ~ |1 -z\ 2 ^~ P + 2 vc- 2pcos0). (6.8) 

The last term is minimal for = 0, where it equals 

(v^-p)(p + v^ + 2) >0, (6.9) 

for p G (0, y/c\. 7 P is a steep descent path for f because the value zero 
is attained only for p = sfc and, in that case, only at one point, = 0. 
However, close to the critical point it is not optimal, because the steepest 
descent path leaves z c at angle ±7r/3 (there are rays where lm(z — z c ) 3 = 0). 
By symmetry, we need to consider only x > 0. 

dRe(f (z = -^Tc + e-^ 3 x)) = x 2 Q(x) 

dx 1 2; | 2 1 1 — z 1 2 ' 

with Q(x) = y/c{l+y/c)— x(l+x)/2— y/cx. Q(0) > 0, and the computation of 
the (at most) two zeros of Q(x) shows that none are in the interval [0, \/c/2]. 
Thus 710c is also a steep descent path for f . Since this is the steepest descent 
path for /o around the critical point, we choose as path r in I n ,t{x) the one 
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formed by 7i oc close to the critical point, until it intersects 7 p _^^ ,and 

then we follow 7 /^jg- 

The Taylor expansions near the critical point z c = —\fc of the functions 
f k are given by 

fo(z) = /o(-v^) + ^o(^+v^) 3 + 0(^ + v^) 4 ), «o= rn \ rv 

y/C{L + y/C) 

f2(z) = f 2 (-V^) + K 2 (z + v ^) + 0((z + v rc) 2 ), K2 = --^, 

1 + V c 

f 3 (z) = -ln(-Vt) + 0(z + yty. (6.11) 
Thus in our case we have 



a = «o = l/(Vc(l + Vc)), 6 = 0, c=-s/(l + Vc), e /s(Zc) = -1/Vc. 

(6.12) 

This, together with the relation 

exp(t/ (z c ) + t 2 "h{z c ) + t^f 2 {z c )) = (±±^¥^ (6.13) 
proves (16. ip . □ 
Lemma 6.2. Fzx £ > and consider the scaling of Lemma \6.1\ Then 



< const e" s , (6.14) 

uniformly for s > —£, where const is a constant independent of t. 

Proof of Lemma \6. 6 A For any finite £, the bound for s G [—£, £] is a conse- 
quence of Lemma 16.11 The value of I can be chosen large but independent 
of t. The strategy for s > i is just a small modification of the computation 
made in Lemma 16.11 and was already used for example in Proposition 17 
of [10] and in Proposition 5.3 of [7]. Let us explain it. 

In Lemma |6~T1 we have seen that 7 P is steep descent path for / , for any 
p E (0, Vc]. Set s = (s+£+£)t- 2 / 3 > £t~ 2/3 > and/ (z) = f (z)+s\n(l-z). 
For any s > 0, 7 P is also a steep descent path for fo(z). However, for s > 
there are two real critical points for / , say at zf with > \%c\- For s 



small, we have at lowest order in s, zf ~ — v^-F Vs ^/k 2 /ko, with k,q and k 2 
given in (16. lip . To get the best bound we should pass through z~ . However, 
this precision is not needed to get exponential bound and we can choose the 
integration path passing through 

-P=l ~^r + ? K2 i K °l / /2 (6-15) 
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With this choice, for e small enough, we have — y/c < —p < z~ and in 
particular, for small s, —p is very close to the position of the critical point. 
As in Lemma 16.11 we use the fact that 7 P is steep descent to control the 
contribution away from \z + p\ < 5, while the contribution close to z — — p 
is controlled by the Taylor expansion of Re(/o(z)), leading to a Gaussian 
bound. By choosing I large enough, all the terms coming from Ke(fk(z)), 
k — 1, 2, 3 are dominated by the leading term of Ke(f (z)). The final result 
is that 

< const Q(p), (6.16) 



e~V a (l + a/c) x 



with 



Q(p) =exp(Re(t(/o(-p)-/o(^))+t 2/3 (/i(-p)-/i(^))+t 1/3 (/2(-p)-/ 2 (^))). 

(6.17) 

Q(p) is decreasing for — p from z c to z c and — p — z c is at most of order 
sfe. Thus, we can easily bound Q(p) by using Taylor expansions. Simple 
computations lead to the desired exponential bound. □ 

To get the needed bound on q n (see f l4.10l) -( pLT4l) ) around the edge, we use 
the bound of Lemma [6.21 on J n 4 which has still to be multiplied by B n> t(x) 
given in (14. 14f) . 

Lemma 6.3. Let n = ct, x = (1 + \/c) 2 t + st 1 ^ 3 and fix i > 0. Then 

\q n {x,t)\ < const r 1/3 e _s , (6.18) 
for any s > — t, and const is a t-independent constant. 
Proof of Lemma \6.3[ This result follows from Lemma 16.21 if 

e -Vct(l i AT) a: 

B n , t (x)= B n>t {x) 1 ; ; < const. (6.19) 

(-Vc) n 

For the factorials we use Stirling formula, namely 

„! = V2^(-) e>», r _^_ </„<_. (6.20) 

We obtain 

B ct , t ((l + V~c) 2 t) = ( (1+ °^ 2 ) V4 (1 + 0(1/*))- (6-21) 
For x = £i, £ G [(1 + v^c) 2 , oo), we compute 



^((i + v^yt) 



( (1+ ^ )2 y /4 (l + 0(l/t))e^, (6.22) 
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with 

h(C) = §£(1 - HO + 2 ln(l + VS)) - |(1 + v^) 2 - (6.23) 

Since h'(£) = at £ = (l + v^) 2 and = -1/(20 < 0, we have e th{ - x) < 1. 

□ 



6.2 Asymptotics in the bulk 

In this section we derive a precise expansion of I n .t{x) for x = At, with 
A G ((1 — \/c) 2 1 (1 + \fc) 2 )- For any fixed c > 0, set n = ct and x = At. Then 



I n ,t( x ) 



2ni 



dw 



r 



w 



(6.24) 



sec 



for the definition of G. Recall a few results from Section [3J For 
A G ((1 — y^) 2 ' (1 + "\/c) 2 )) 9 has two complex conjugate critical points, w c 
and w c , with w c = f2(A,c, 1). In particular, \w c \ = y/c, |1 — w c \ = vA, and 
|g"(w7 c )| = ^^a/4c — (1 + c — A) 2 . When (r/, ^, r) = (c, A,l), we denote by 
7r c the angle and by 7Ta the angle n u . Then 



t-./ / s\ 1 -H c — A c , . . A,.,. 

RefoM) = 2 2 ln(c)+ 2 ln(A) ' 

lm(g(w c )) = lm(w c ) - Att c - ctt a , 

arg(a(w; c )) = ~ + tt c - tt a . 



(6.25) 



Lemma 6.4. Set a = a(c, A) = lm(g(w c )) and f3 = /3(c, A) = — \{k c + 7r A + 
7r/2). T/ien, as t — > oo ; 



/ct,*(At) 



e tRe(g(ui c )) 



7T W c 



cos(to + /?) + C(r 1/2 ) 



(6.26) 



For any sq > 0, the errors are uniform for X G [(1 — \/c) 2 + £o> (l + \/c) 2— £ o]- 
Remark 6.5. In fact, we prove the bound of A6'.26}) with error term 

0{t- 1 ' 2 ) + 0(vV(«; c )|te- const l9 " (Wc)|52 *) (6.27) 



for some < 5 <C |o"(w c )|. In Lemma \6. 7| we will have to be careful with the 
second term of the bound, since g"{w c ) goes to zero at the edge. 



82 



Proof of Lemma \6.4\ The critical points of g, the points such that g'(w) = 0, 
are w c and its complex conjugate w c . Close to w c the Taylor expansion of g 
has a first relevant term which is quadratic, 

g{w) = g(w c ) + \g"{w c ){w - w c f + 0{{w - w c f). (6.28) 

Now we construct the steep descent path used in the asymptotics. By 
symmetry we consider only Im(u') > 0, the path for Im(w) < will be the 
complex conjugate image of the first one. Let 7 P = {w = pe lct> ,(f) G [0,7r]}, 
then 



^-(R e (g(w = pe>*))) = psm(4,) 



X -1 



1 — w\ 



(6.29) 



This is positive if |1 — w\ < \f\, and negative otherwise. 

Locally, consider the path 7i oc = {w = w c + 9x, x G [—5, 5}}. Then 



g(w) = g(w c ) + \g"(w c )e z x 2 + 0{x A ), (6.30) 

where we choose 

= ex P - \ arg(/(w c ))^ = exp + l J^_lA^j . ( 6 . 3 i) 

For — 5 < x < 0, the path 7i oc is closer to 1 than \^X, while for < x < 5 
the path 7i oc is farther from 1 than vX This is the case since our 7i oc has 
an angle between 7r/4 and 37r/4 to the tangent to the circle |1 — w\ = vX 

So, the steep descent path used is the following: we extend 7i oc by adding 
two circular arcs of type j p , for adequate p, which connect to the real axis; 
finally we add the complex conjugate image, see Figure IB~2l too. In this way, 
we have a steep descent path. Thus, 

I nt ( x ) = e *Re( 9 K)) ( e -^) +2 Re(— [ ^ e t9(w) ) (6.32) 

V 2 W 7loc w J 

with /i ~ \g"(w c )\d~ 2 , as soon as \g"(w c )\ > 0, i.e., as soon as the second order 
term dominates all higher order terms in the Taylor expansion. 
The second term of f 1 6 . 3 2 [) is given by 



- I — e«M = — f dx-e^e-^'^e^Oix) 

r-<5 



— — [ dxe tg{^) e -lt\ g "( Wc )^ + Ei / 6>33 x 

7T1 Wc J -5 
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Figure 6.2: Illustration of the steep descent path. 



where 



E l = — — I dxe t9 ^e-\ tw ' {w ^ x2 e ^0{tx\x). (6.34) 

27T1 W c J_ s 

Here we used \e x — 1| < |x|e' z '. Changing variable y = xy/i, we get that 

\Er\ < conste ffle ^»- / dye~\ 9 "^^ 2 0(y } y 3 /Vt)e°^^ 

t JsVt 



e tRe(g(w c )) 



< const — . (6.35) 

for 5 small enough, i.e., for < 5 <C \g"(w c )\. In this small neighborhood, 
the quadratic term controls the higher order ones. The final step is to extend 
the integral on the rest of r.h.s. of (I6.33P from ±5 to ±oo. This can be made 
up to an error e iRe ^^ c ^O(e~ Mi ) as above. 

Resuming we have (counting the contribution from both critical points) 



J f _ e tRe(g{w c )) 



0{e->*) + 0(1/(«VI7WI) 

+2R e ( —A I dxe t 9 M e -W'^\A . (6.36) 
\2ni wJ R J 1 ; 

The error terms are the one indicated in (I6.26p . and the Gaussian integral 
for the last term gives 

2f " U " 1 " " Re f-i^ e i^W»c))] . ( 6 .37) 



^/27it\w c \ 2 \g"(w c )\ V w 
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We then set (3 = arg(-i#/w c ) = -vr/4 - (vr c + tt a )/2, so that -i0^ = e i/3 . 
For A in a compact subset of ((1 — \/c) 2 , (1 + y/c) 2 ), \g"(w c )\ is uniformly 
bounded away from zero and infinity. Thus the Lemma is proven. □ 

The bound of Lemma 16.41 can be easily extended until position away of 
order 0{t 2 / 3 ) from the upper edge. 

Lemma 6.6. Set a = Im(gr(itJ c )) and (3 = — |(vr c + n\ + 7r/2). Then, for 
A G [(1 — y/c) 2 + s , (1 + \[c) 2 — t^ 1 / 3 ], for any fixed e > 0, we have the 
uniform estimate 



e tRe(g(w c )) 



V\9"M\t 



cos(to + 0) + 0(r 1/2 ) + O (Vie 



-const t 1 / 3 

7T\W C \ Z V 

(6.38) 

Proof of Lemma \6.6l The analysis of Lemma 16.41 can be made also for this 
case, with only minor differences. Indeed, for (1 + \fc) 2 — A ~ t^ 1 ^ 3 , we have 
\g"(w c )\ ~ t -1 / 6 and this time we choose S going to zero as t — > oo, setting 
5 = t -1 / 4 . With this choice, (16.321) and (16.351) are still valid because at the 
border of integration the quadratic term dominates the cubic one. Indeed, 
with y = 5Vi = t 1/4 , it holds y 3 /Vi ~ t 1 / 4 < t 1 ' 3 ~ \g"{w c )\y 2 . Also, the 
error term coming from steep descent in (I6.32|) will vanish as t — > oo but 
slower than before, with [it ~ t 1 ^ 3 . □ 

The results of Lemma 16.41 and Lemma l6~o1 imply the following asymptotics 
for the functions g n 's. 

Lemma 6.7. Set a = lm(g(w c )) and (3 = — \{n c + 7r A + tt/2) and fix any 
e > 0. Then, uniformly in A G [(1 — y/c) 2 + e , (1 + y/c) 2 — t~ x ^\, we have 



1 r 1 / 2 

q ct (Xt,t) = 



7 N c _ (!+ c - A ) 2 



cos(ta + /3) + 0(t~ 



-l/2\ 



(6.39) 



Proof of Lemma \6. 1\ We just have to compute the pref actor 
B ct! t(\t)e tRe ( g ( Wc ^ . We have (I6.25P and applying Stirling formula for the fac- 
torials in B ctit (\t) we get that 

B ct ^\t)e me{9{w ^ = (A/c) 1/4 (l + 0(l/t)). (6.40) 

□ 

Now we need to fill the gap between the bulk and the edge. In this region 
we do not need precise asymptotics, just a bound. Approaching the upper 
edge, g"(w c ) goes to zero, but then everything can be controlled by the cubic 
term, because \g"'(w c )\ ^ at the edges. 
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Lemma 6.8. For sq > fixed but small enough, and £ > large enough, we 
have the bound 

£-1/2 

\q c t(Xt,t)\ < const (6.41) 

4/~ (1+C-A)2 



C 4 



uniformly for A G [(1 + v/c) 2 - e , (1 + v^) 2 - £*~ 2/3 ]. 

Proof of Lemma \6.8l Consider the eo-region close to the upper edge with 
eo > small enough. We can compute explicitly the direction 8, see (16.311) . 
It is a continuous function of A and, as A | (1 + v^c) 2 ? $ T e l5?r / 4 (because 
n x t 7t and 7r c | 0). We need just a bound, so we choose 6 = e l57r//4 and set 
the local path as 



Tic 



{ w = Wc + e iW4 x, x G [-5, Im(w c ) V2]}. (6.42) 



The path 7i oc reaches at x = Im(«; c )\/2 the imaginary axis and this is the 
reason for the upper edge of 7i oc . We have 

g{w) = g{w c ) + \g"{w c ) (w - w c ) 2 + \g'"{w c ) (w - w c ) 3 + 0{{w- w c ) 4 ). (6.43) 

In a 5-neighborhood of w c , along the direction 8 chosen, 

Re(lg"(w c )(w - w c f) = -\\g" {w c )\x 2 {l + 0(e o )) (6.44) 

and 

Re(±g'"(w c )(w - w c f) = l\g"'(w c )\x 3 /V2(l + 0{^)). (6.45) 

Therefore, for 6q small enough, the quadratic term helps the convergence. 
For x < 0, the cubic term helps the convergence, while for x G [0, lm(w c )V2] 
we will need to control it by the quadratic term. Thus, 



I n t ( x ) = e *Re0K-c)) o(e -M*) + 2Re (— [ ^aM) 

V27ri^ oc w J 



(6.46) 



with fi ~ \g"'(w c )\5 3 , where g"'(w c ) 2/^(1 + y/c) as A -> (1 + v^) 2 . 

Consider then the contribution coming from the integral over 7i oc . We 
have 



dw c tp(w) 



w 

floe 



< 



e ffle(j(w c )) />Im(w c )V2 



Wr 



dxexp (-\t\g"{w c )\x 2 ) (6.47) 



x exp (M\cf"{w c )\x z /y/2 + O(xH)) (1 + 0(x)), 
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the last l/v2 coming from Re(e 17r / 4 ) = 1/V2. A simple verification gives 

- \t\g"{w c )\x 2 + \t\g"'{w c )\x 3 /y/2 < -\t\g" \w c )\x\ < x < lm(w c )V2. 

(6.48) 

So, for x E [0, lm(w c )y/2], the quadratic term is still dominating higher order 
terms, including the cubic one (the quartic term can be bounded by replacing 
\ by | in the above estimate). 

On the other hand, for —5 < x < 0, we have that the cubic term is 
negative and dominates all higher order terms. More precisely, for 5 small 
enough, 



exp (^t\g"'(w c )\x 3 /y/2 + 0(xH)^j < exp (±t\g"'(w c )\x 3 ) < 1, (6.49) 

in the region x G [—5, 0]. 

Using (I6.48P for positive x and (I6.49P for negative x, we get 

dxexp {-\t\g"(w c )\x 2 ) 

< m ns tP fflc (gK)) 1 - (6.50) 

Replacing the value of \g"(w c )\ into this expression ends the proof. □ 



6.3 Asymptotic of the kernel 

In this section we obtain the precise asymptotics of the extended kernel in 
the bulk first and a bound to control the behavior starting from the upper 
edge. Here we use several notations introduced in the Section [3j 

As usual, it is convenient to conjugate the kernel before taking the limit. 
For the upper edge (cf. Lemma 16.21) set 

W itU = exp ^-Vn~Ti + Xi ln(l + ^/n~/U) - ln(— \/nJti) - tij , (6.51) 

and, in the bulk (see Lemma 16.41) set 

W ijh = exp (|(ti + rii - Xi) - \ni\n{ni/ti) + ~Xj ln^j/ij) - U) . (6.52) 

Then, define the conjugation as 

w = [Wi^ for (y/U- ^/Wi) 2 < Xi < (v^+ ^TiT) 2 , , . 

1 \W i>u , for Xi > (V^+ y^") 2 . l ' ; 

Remark that Wi is continuous. Moreover, \Wi jU — = 0(L _1 / 3 ) for 

\xi — (y/U + v^") 2 ! = OiL^ 1 / 3 ). Therefore in such a neighborhood it is ac- 
tually irrelevant which formula to use. 
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Proposition 6.9. Let us consider two triples (xi,7ii,ti) and (x 2 ,n 2 ,t 2 ) pa- 
rameterized by 

%i = [uiL], Hi = [rjiL], U = TiL. (6.54) 

Assume that they are in the bulk of the system, namely, that exists 6q > 
such that 

(Vn - Vv^) 2 + eo<Vi<(Vn + Vvl) 2 - L ~ 1/3 - (6-55) 

Denote z c = fi(i/i, 771, 7i), w c = VL{y 2l rj 2l T 2 ) , and assume that these points 
are not too close: \z c — w c \ > L -1 / 16 . Then, the asymptotic expansion 

{W 1 /W 2 )K(x 1 ,n 1 ,t 1 ;x 2 ,n 2 ,t 2 ) 1 



2ixL^\G"{w c )\\G"{z c )\\l 



z, 



X 



e iLIm(G(u> c ))+i/3 2 e iLIm(G(w c ))+i/3 2 

+ ~ ~ —\Tl™<a<r-^-\R, (6.56) 



Wc — Zc giLIm(G(z c ))+i^i yj^ _ 2 C g— iiIm(G(z c ))— i/3i 

I e -iLlm(G(w c ))-i/3 2 1 g-iiIm(G(«j c ))-i/3 2 ^ g 

+ ^TT^ e iLIm(G(^ c ))+i / 3i + ^ c _ ^ e -iLIm(G(^ c ))-i/3i + > 

holds, with the error uniform in L for L > L Q ^> 1. TTie phases Pi and f3 2 
are given by 

A--T~"2 2~' ^ 2 "T + 1 2~" (6 ' 57) 

Proof of Proposition \6.9l The analysis relies on the double integral represen- 
tation (14.21) of the kernel. The analysis for the cases (ni,ti) 7^ (n 2 ,t 2 ) and 
(rii,ti) -< (n 2 ,t 2 ) are very similar. Let us explain the first case, correspond- 
ing to 771 > 772, or Ti < r 2 , or (771,71) = (772,72). The asymptotics employs 
several ingredients already used in Lemma 16.41 and Lemma 16.61 Thus, we 
introduce the notations 

Ci = r)i/Ti 7%i= [ciU], Xi = Ui/n Xi = [XiU]. (6.58) 

The conjugation factor e* 1- * 2 in the kernel representation (14.21) will not appear 
in the following computations, since it appears automatically in the factors 
Wi/W 2 . Thus, we have to analyze 

W?L iw L d ^" Ml " w (!-,)(<»-,) (6 ' 59) 

with gi(w) = w + Xi ln(l — w) — Ci ln(w) = G(w\Xi, Ci, 1), i = 1, 2. 
The critical points of g 2 (w) and gi(z) are given by 

w c = Q(X 2 ,c 2 , 1) = n(i/ 2 ,T72,r 2 ), z c = fi(Ai,ci,l) = fi(vi, 771,71). (6.60) 
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Figure 6.3: Illustration of the steep descent paths. 




The integrals over w are, up to the factor w/(z — w), as in Lemma 16.41 
Therefore, the steep descent path T is chosen as in Lemma 16.41 and the 
steep descent path Ti is chosen in a similar way. We illustrate these paths 
if the critical point is ( c , see Figure IBTBl In particular, |to c | = ^/7/2/t 2 and 
\ z c\ — \A?i/ r 2- m our case, we have \w c \ < \z c \ and \w c — z c \ > L -1 / 16 . The 
steep descent paths described above actually intersect. Therefore, we have to 
correct (16.591) by subtracting the residue at z = w, as indicated in Figure [6T4l 
We call "main term" the contribution of the integral with To and Ti crossing, 
while we call "residual term" the contribution of the residue. 

Notice that the integral with the paths r and Ti crossing is integrable 
in the usual sense, because the divergence term l/(w — z) is integrable. The 
contribution of the main term is the following. 

Both integrals can be divided as the part in H and its complex conjugate. 
Therefore, in the final expression we get the sum of four terms. Now, we 
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restrict our attention to the integral over the path T and Ti on H. The 
analysis of the integral over T is the same as in Lemma 16.41 except for the 
missing l/w c factor and that instead of 2Re(- • • ) we just have (• • • ) in (16.361) . 
The integral over I\ is similar. 

This time we choose the cutoff for the evaluation of the term with the 
steep descent path equal to S = L~ 1//4 . There are two reasons. The first 
one is that we want to get the expansion valid also for Vi up to L -1 / 3 away 
from the upper edge, compare with Lemma 16.61 The second reason is that 
we have the extra factor l/(w — z). The contributions of the steep descent 
path do not create problems, since the factor is integrable in the usual sense 
(just need a bound). However, with our choice of 5, in the contribution of 
the 5-neighborhoods of z c and w c we have 



l/(w -z) = l/(w c - z c ) + 0(5/\w c - z c \ 2 ) = l/(w e - z B ) + 0(L'^), (6.61) 

with 5 = L" 1 / 4 and \w c — z c \ > L -1 / 16 . In the end, the contribution of the 
main term is given by 



D t 2 Re(g 2 («'c)) 



w c - z c ^2itt 2 \g'i(w c )\ 



^ lm{gM) 9 2 (w c ) + O (L- 1/2 ) 



-tiRe(gi(z c )) 



X 



+ 



^2Tit l \g'{(z c )\ \l-z c 



e t2~Re(g 2 (w c )) 



-tiRe( S i(z c )) 



y/27rt 2 \g%(w c )\ y^\gj(zj\ |1 - z c 



rO{L 



-1/8N 



(6.62) 



The term e 17Tci is the phase of 1/(1 — z c ), while 9i are the directions of the 
steepest descent paths at the critical points. Explicitly, 

9i(z c ) = exp(i(7r - \ axg(g"(z c )))) = exp(i37r/4 + i(vr Al - 7r Cl )/2), (6.63) 
9 2 {w c ) = exp(i(| - | axg(5f"(«J c )))) = exp(i57r/4 + i(7r A2 - vr C2 )/2). 

Putting together the four terms (two times two critical points) we get the 
complete contribution of the main term as 



3 t 2 Re(g2(wc))-tiRe(3i(2 c )) 



2irWtit 2 \l 



Zr. 



\g'2>o)U{z c )\ 



3 it2lm(s2(Wc))+i/32 



+ 
+ 



+ 



e it2lm(g2(w c ))+if3 2 



Wc 



z r e 



-it2lm(32(«>c))-i/32 



-itilm(gi(z c ))-i/3i 
3 — it2lm(g2 (w c )) — i/?2 



yj — z gitilm(gi(« c ))+ift 



+ 



W r 



z r . e 



-itilm(gi(z c ))-i/3i 



(6.64) 
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with Pi = — arg(0i(z c )) — n Cl and p 2 = axg(9 2 (w c )). Finally, we replace 
gi(w)ti = G(w\v il rj il Ti)L, tt Ai = ir u ., n Ci = and 

e t 2 Re(g 2 (w c ))-t 1 Re(g 1 (z c )) e t 1 -t 2 _ W 2 /Wi (6.65) 

to get (1535)1 . 

The final step is to estimate the contribution of the residual term (the 
last case of Figure 16.41) . It is given by 

1 e ( T 2- T i) Lz e (Vi-V2)Lln(z) 

2^i J c dZ (l-z)^-^ • ( 6 ' 66 ) 

where £ and £ are the two intersection points of the steep descent path r 
and Ti. Since t 2 —Ti > 0, r\i—r\ 2 > 0, and |1 — z| = const along the piece of 
r x inside r , we have Ke(z) < Re(C) and Re(ln(z)) < Re(ln(£)). Therefore, 

| (16.661) | < e *2R- e (92(C))-*iRe(si(C)) < e *2Re(52(ui c ))-tiRe(ffi(* c ))^^ e -/iiti e -/i2*2^ 

(6.67) 

for some positive /ii, fx 2 and at least one larger than L -1 / 8 . This follows from 
the fact that either one (or both) critical points are away of order L -1 / 16 
from (, and ( lies on the steep descent paths of g 2 {w) and —gi(z), with local 
quadratic behavior. □ 

While doing time integration we will also need the following corollary. 

Corollary 6.10. Consider the same setting of Proposition [679\ Then 

(a) the formula for K(xi, nx, ti,x 2 + l,n 2 ,t 2 ) is the same as k6. 56]) but with 
an extra factor (1 — w c ), resp. (1 — w c ), to the terms with e 1 ^ 2 , resp. e~ 1/32 . 

(b) the formula for K(xi,rii — 1, t±; x 2 , n 2 , t 2 ) is the same as $6.56]) but with 
an extra factor z" 1 , resp. (5 C ) _1 ; to the terms with e~ 1/Bl , resp. e 1/Sl . 

Proof of Corollaru \6.10\ The proof is almost identical to the one of Proposi- 
tion 16.91 The only difference is that in (16.591) we have for (a) an extra term 
(1 — w) and for (b) an extra 1/z. □ 

At this point we have all the needed estimates in the bulk. However, 
since our system develops facets, we need to have control at the upper edge. 
We will just need some bounds and, since the integrals are the same as in 
Section [6J], apart from the factor l/(w c — z c ), which we will assume bounded 
away from zero. 

Proposition 6.11. Consider the setting of Proposition l679\ but with one or 
both of the Vi close to the upper edge, 

(Vn + v 7 ^) 2 - L ~ 1/3 <"*<(Vn + V^u) 2 - ^~ 2/3 - (6.68) 
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Moreover, assume that \z c — w c \ is bounded away from zero uniformly in L. 
Then, there exists £ large enough, such that 

const 

(W 1 /W 2 )\K(x 1 ,n 1 ,t 1 ;x 2 ,n 2 ,t 2 )\ < . _ (6.69) 



L Ui=i \ vm - Hn + Vi 



uniformly in L for L > L 1 . 

Proof of Proposition \6AT\ The proof follows the same argument as Lemma I6~B1 
for the variables which are close to the edge. For the one which is away from 
the edges, it is a consequence of the analysis Proposition 16.91 □ 

When one or both positions are at the edge, we need a different bound. 

Proposition 6.12. Consider the setting of Proposition f^TPL but now with v 2 
at the edge or in the facet, i.e., 



for any fixed i. Assume \z c + \Jrj 2 ~]t 2 \ is bounded away from zero uniformly 
in L. Then, 

coiisti 

(Wi/W 2 )| J R:(xi,n 1 ,ti;x2 > n 2 ,t2)| < 1 (6.71) 

v^yvin - \{n +r]i - z/i) 2 

1 ( x 2 - {y/T^ + ^m) 2 L 



- L i/3 ex P^ (r 2 L)V3 
uniformly in L for L > L 1 . 

Proof of Proposition \6.1S[ The proof is obtained along the same lines as 
Lemmas 16.11 and 16.21 With respect to those cases, the integral has however 
an extra factor l/(w — z). Since we need just a bound, it can simply be 
replaced by l/(w c — z c ) as follows. In Lemma EH] w c is replaced by —^/rj 2 /T 2 , 
while in Lemma 16.21 we need to replace w c by p as given in fl6.15p . Notice 



that in the last case we can take \w c + y^/^l as small as desired. The 
assumption \z c + \rrj 2 /r 2 \ > uniformly in L ensures then that l/(w c — z c ) 
remains bounded as L —>■ oo. □ 

The last case to consider is when both v\ and v 2 are at the upper edge. 

Proposition 6.13. Consider the setting of Proposition tfTQ but now with V\ 
and v 2 at the edge or in the facet, i.e., with 

Vi>{^n + Vni) 2 -ZL- 2 '\ z = l,2, (6.72) 
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for any fixed £. Assume that | \Jrj 2 ]T 2 ~~ \ph]~T\\ is bounded away from zero 
uniformly in L. Then, 

(W 1 /W 2 )\K(x 1 ,n 1 ,t 1 ;x 2 ,n 2 ,t 2 )\ < const (6.73) 

LV^ XP \ (r 2 L)V3 jexPy (riL) i/3 

uniformly in L for L > L ^> 1 . 



Proof of Proposition \6.13[ The proof is like Proposition 16. 121 We will have 
\ w c + \A? 2 / r2 l anc ^ l z c + ^/vi/ T i\ as small as desired. The assumption 
I \A?2/' r 2 — y^hh\\ > uniformly in L allows us to easily bound uniformly 
in L the term l/(w c — z c ). □ 



A Determinantal structure of the correlation 
functions 

Let £i, . . . , Xjv be finite sets and c(l), . . . , c(N) be arbitrary nonnegative 
integers. Consider the set 

X = (Xi U • • • U Xi) U • • • U (Xat U • • • U Xat) (A.l) 

with c(n) + 1 copies of each X n . We want to consider a particular form of 
weight W(X) for any subset Id, which turns out to have determinantal 
correlations. 

To define the weight we need a bit of notations. Let 

<p n ( ■ , • ) : X n _i x X„ C, n = 2, . . . , N, 

n (virt, • ) : X n -»• C, n = 1, . . . , N, (A.2) 
:Xiv^C, j = 0,...,iV-l, 

be arbitrary functions on the corresponding sets. Here the symbol virt stands 
for a "virtual" variable, which is convenient to introduce for notational pur- 
poses. In applications virt can sometimes be replaced by +oo or — oo. The 
4> n represents the transitions from X n _i to X n . 
Also, let 

t N < ... <+ N - t 1 *- 1 < ■■■ < z-^- 1 - t N ~ 2 < ... <t 2 - t 1 < ... < t 1 

£(3 — — l c(N) — l ^ ^ l c(N-l) ~ l — — l c(2) — l ^ - - - ^ t c (l) 

(A.3) 

be real numbers. In applications, these numbers refer to time moments. 
Finally, let 

^*S.*a_i( ' ' ' ) : %™ x %n —> C ' 71=1, ...,iV, a = 1, . . . ,c(n), (A.4) 
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be arbitrary functions. The %n jt n represents the transition between two 
copies of X n associated to "times" t n a _ x and t n a . 

Then, to any subset X C X assign its weight W(X) as follows. W(X) 
is zero unless X has exactly n points in each copy of X n , n = 1, . . . , N. In 
the latter case, denote the points of X in the mth copy of X n by a;]J(t^J, 
k — 1, . . . , n, m — 0, . . . , c(n). Thus, 

X = K(C) I k = 1, . . . , n; m = 0, . . . , c(n); ra = 1, . . . , N}. (A.5) 

Set 



n (A ) = H 



det W^O, *?(**»)))] 



KfcKn 



71=1 

c(n) 

X 



ndet [^.^^(o.xrcci))]^ 

a=l 



det W-,^^))]^,^, 
(A.6) 



where • ) = virt for all n — 1, . . . , N. 

In what follows we assume that the partition function of our weights does 
not vanish: 

Z := W ( X ) ^ °- ( A - 7 ) 

xcx 

Under this assumption, the normalized weights W(X) = W(X)/Z define a 
(generally speaking, complex valued) measure on 2 X of total mass 1. One 
can say that we have a (complex valued) random point process on X, and its 
correlation functions are defined accordingly, see e.g. [15]. We are interested 
in computing these correlation functions. 

Let us introduce the compact notation for the convolution of several tran- 
sitions. For any n — 1, . . . , N and two time moments t™ > t% we define 

Tin fn = Tin fn >K Tin fn >[:•••* Tin fn sj~n _ /J" (A. 8) 

where we use the notation (/ * g)(x,y) := J2 Z f( x ' z )d(, z i v)- F° r an y time 
moments > with 7^ (^2,^2), we denote the convolution over 

,n2 \ 

all the transitions between them by (fP a i-> ta 2>\ 

fx 71 ! j- n 2 \ 11 

(t ai ,t O2 ) = * 0ni+1 * T n 1+ 1 * . . . * * ^ 2))t:| . (A.9) 

If there are no such transitions, i.e., if < t™'* or (ai,ni) = (a 2 ,n 2 ), we set 

0(*ai,*a|) — Q 
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Furthermore, define the matrix M = ||M^|L , =1 by 

M Kl = (fa * T k * ■ • ■ * <j> N * T N * (virt) (A. 10) 

and the vector 

*:1=0 ( ^ } **£U i = i,...,n. (A.n) 

The following statement describing the correlation kernel is a part of Theorem 
4.2 of [7]. 

Theorem A.l. Assume that the matrix M is invertible. Then 
Z = detM 0, and the (complex valued) random point process on X de- 
fined by the weights W(X) is determinantal. Its correlation kernel can be 
written in the form 

K(Q,x 1 ;Q,x 2 ) = -^^i\ Xl ,x 2 ) 

N n 2 

+ E E <l'S MiM- 1 )^ * ( ^)^ 2 ))(vfrt, x 2 ). 

k=l 1=1 

(A.12) 

The proof of Theorem IA. II given in [7] is based on the algebraic formalism 
of [15]. Another proof can be found in Section 4.4 of [26]. Although we stated 
Theorem I A . 1 1 for the case when all sets X n are finite, one easily extends it to 
a more general setting. Indeed, the determinantal formula for the correlation 
functions is an algebraic identity, and the limit transition to the case when 
X n 's are allowed to be countably infinite is immediate, under the assumption 
that all the sums needed to define the *-operations above are absolutely 
convergent. Another easy extension (which we do not need in this paper) is 
the case when the spaces Xj become continuous, and the sums have to be 
replaced by the corresponding integrals over these spaces. 
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